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Abstract 


A detailed study of the spectroscopy and excitation dynamics of the trivalent 
lanthanides Tm 3+ and Ho 3+ in Yttrium Lithium Fluoride, LiYF4 (YLF), has been done. 
YLF is a very versatile laser host that has been used to produce laser action at many 
different wavelengths when doped with trivalent lanthanides. Since the early 1970’s YLF 
has been the subject of many studies, the main goal of which has been to produce long 
wavelength lasers in the eye safe 2|im region. This study concentrates on a presentation 
and analysis of the spectroscopic features, and the temporal evolution of excitation 
energy in YLF crystals doped with Tm 3+ and Ho 3+ ions. 

Absorption spectroscopy is studied to identify wavelength regions where energy 
can be absorbed in Holmium YLF and Thulium YLF, and to determine their respective 
absorption cross sections. These measurements are applied in the Judd-Ofelt theory to 
determine radiative transition rates of spontaneous emission. Luminescence spectroscopy 
is studied under cw diode laser excitation at 785nm. The effect of dopant ion 
concentration and excitation power on the observed luminescence are considered in these 
measurements. An analysis of these measurements have been used to determine channels 
of energy transfer between Tm 3+ and Ho 3+ ions (cross relaxation, upconversion, and 
resonant energy transfer). The temporal response of Tm and Ho in singly and co-doped 
YLF to pulsed laser excitation with a Ti:Al203 laser and a CoMgF2 laser tuned to various 
wavelengths have also been studied. The energy transfer mechanisms of cross relaxation, 
upconversion, and resonant energy transfer between Tm 3+ and Ho 3+ ions have been 
modeled, and the model parameters extracted by a fitting procedure to the measured 
temporal response curves. Rate equation approaches to modeling are presented that result 
in predictions of rate constants for energy transfer processes, as well as more 
conventional approaches to modeling such as the Forster-Dexter models, which give the 
interaction strengths in terms of microscopic interaction parameters. 
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CHAPTER 1 

Introduction 


Every story involves one or more archetypes. To make 
a good story a single archetype is usually enough. 

Umberto Eco 
Travels in Hyperreality 

1.1 Overview 

A notable feature of the rare earth (RE) ions when doped into crystals is the large 
number of sharp lines in the absorption and emission spectra. The wavelengths of these 
lines range from the ultraviolet to the infrared. Solid state crystals, which are insulators, 
have a band gap that corresponds to the energy of an ultraviolet photon, so that visible or 
infrared luminescence is not expected from the pure (undoped) crystal. The bulk material 
in which the RE ions are imbedded will have a perturbative effect on the ions that serves 
to shift and split the free ion energy levels. These perturbed electronic states can also 
interact with the phonon states of the crystal lattice to produce various effects. Despite 
this, the situation is greatly simplified since the band structure of the solid does not play 
a part in the absorption and emission characteristics of the RE ion. 

An emitting solid is, nevertheless, a complicated system consisting of a large 
number of ions, each of which is capable of interacting with its neighbors. Once the ions 
are excited with a suitable source, the effects of these interactions becomes manifest. 
Generally, in the case of RE ions a diverse assortment of luminescent phenomena are 
observed. The challenge is to differentiate between the various luminescent behaviors that 
are observed, to identify the channels through which energy is transferred, and to 
determine the nature of the interaction that promotes transfer. 



This situation is complicated enough in materials that contain only a single ion 
species, so it is no surprise that the complexity is increased when ions of different species 
are introduced. This increased complexity, however, has a beneficial aspect in that it can 
be used to produce luminescence in an ion of one species by exciting an ion of another 
species (sensitized luminescence). This, of course, makes possible a diverse range of 
lasing action in RE doped crystals under various excitation schemes. 

This thesis describes a series of experiments designed to probe into the 
spectroscopy and excitation dynamics of Lithium Yttrium Fluoride LiYF 4 (YLF) 
containing Tm 3+ and Ho 3+ ions. The research addresses the full problem of the dynamics 
of excitation energy. This will take into account the spectroscopic details of absorption 
and emission lines in co-doped and singly doped YLF crystals, the channels of their 
excitation and their respective decay behavior. Various experimental parameters have been 
varied, such as concentrations of the dopant ions, wavelength of the excitation source and 
power of the excitation source. This is useful as an aid to differentiate between various 
competing processes which occur in these complex systems, and to facilitate the task of 
determining which of the dopant ions are responsible for which processes, how the ions 
interact with each other, and their relative importance in the overall behavior of the 
response of YLF:Tm,Ho systems to an excitation source. Emphasis is placed on a 
coherent presentation of the experimental data. Careful attention is given towards the 
interpretation of the experimental results and theoretical analysis of the experimental data. 

L2 Historical PacKground 

The early research on YLF based systems concerned 2p.m laser action of 
sensitized Ho 3+ YLF. Remski, et al. [Re69] demonstrated energy transfer assisted pulsed 
lasing action at 77°K of Ho 3+ in YLF :Er,Ho. Chiklis, et al. [Ch71] demonstrated 
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flashlamp pumped high efficiency room temperature laser action at 2.065pm in 
YLF:Er,Tm,Ho. The analysis of the optical spectrum and determination of the energy 
levels were reported for Tm 3+ in YLF by Jens sen, et al. [Je75] and for Ho 3+ in YLF by 
Karayanis, et al. [Ka76], 

The availability of mid-infrared diode lasers in the I980's initiated a resurgence in 
Ho 3+ sensitized lasers because diode pumped operation offers an attractive alternative to 
lamp-pumped schemes and also because the GaAlAs diode emission has a good spectral 
overlap with the absorption of Tm 3+ , making for an efficient, site selective pumping 
scheme. A sampling of studies using diode-laser pumping to achieve 2pm laser action in 
Ho 3+ sensitized YLF include Kintz, et. al. [Ki87], Hemmati [He89], Esterowitz [Es90], 
Budini, et al. [Bu92] and Storm [St93], 

Despite the abundance of studies involving lasing schemes in YLF systems, there 
remains a paucity of spectroscopic studies in the literature. Brenier, et al. [Br89] in an 
extensive paper, studied the excited state dynamics of Tm 3+ ions in YLF and the Tm 3+ -> 
Ho 3+ energy transfer in YLF. In this study they examined the processes of upconversion, 
cross relaxation and energy transfer to explain the observed response of YLF:Tm and 
YLF:Tm, Ho systems to pulsed excitation at various wavelengths. Ozen [Oz91] in her 
Ph.D. thesis examined the energy transfer process between Tm 3+ and Ho 3+ in YLF. This 
study examined the temperature dependence of the decay of Tm 3+ and Ho 3+ in YLF 
under excitation from a tunable dye laser and the spectroscopic features under tungsten 
lamp excitation. 

1-3 Purpose 

Rare earth doped solid state materials have become an important class of solids 
attracting much attention among researchers, as is evidenced by the abundance of studies 
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that can be found in the literature. There are two reasons promoting an incentive to study 
these materials. The first reason is technical and the second rests on scientific grounds. 

From a technical point of view, rare earth materials have been used in a wide 
variety of optical technologies (i.e., lasers, IR quantum counters, television, and lamp 
phosphors.) The research at NASA Langley Research Center has focused on the 
development of reliable, compact solid state lasers in the eye safe 2pm range for remote 
sensing of the earth's atmosphere from space. From a scientific point of view, a basic 
understanding of excitation transport in inhomogeneous and disordered solids has the 
potential for usefulness in many branches of physics. 

As was mentioned earlier in this chapter, a substantial amount of work has been 
done on the lasing characteristics of RE solid state materials. It is the case, however, that 
spectroscopic measurements are somewhat sparse in light of the importance some of the 
solid state lasers have achieved. This has especially been the case regarding Tm 3+ and 
Ho 3+ singly and co-doped YLF systems. This situation seems to be improving. During 
the course of the research for this thesis YLF has seemed to be attracting even more 
attention from a spectroscopic point of view. Most of these studies have concentrated on 
specific topics, the results of which have been presented at scientific conferences in the 
form of slim articles. A comprehensive characterization YLF:Tm,Ho from a spectroscopic 
point of view has not yet been achieved. Many important questions remain unanswered. 
The intent of this thesis is not to answer all of these questions, but to conduct a study 
that covers many topics with the intent of connecting the various aspects in a uniform 
and consistent manner. 

As this thesis will show, especially in the later chapters, YLF:Tm,Ho is an 
extremely complex system. It displays a bewildering plethora of complex and competing 
processes. These processes comprise the topic of energy transfer. Energy transfer 
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manifests itself in various forms (i.e., upconversion, cross relaxation, and direct transfer). 
These various forms of transfer will be explained at such time when they arise in the 
context of the thesis. 

A number of spectroscopic properties are of importance for the modeling and 
understanding of YLF:Tm,Ho lasers. These include radiative decay rates, absorption and 
emission cross sections, energy levels, Boltzmann fractions, and various energy transfer 
rates. This thesis seeks, as one of its goals, to determine specific values for these physical 
parameters whenever possible. This set of values will be useful for those researchers 
attempting to understand and model the optical and temporal behavior of the emission 
from YLF:Tm,Ho. 

1.4 Preview of Contents 

Spectroscopic measurements of the absorption in YLF:Tm and YLF:Ho were 
made, and used in the Judd-Ofelt theory to calculate radiative transition rates between 
sets of energy levels (manifolds). Radiative transition rates were also calculated 
independently by comparison of measured emission cross sections and emission cross 
sections derived from absorption cross sections. 

Luminescence measurements under cw diode excitation have been made. The effect 
of dopant ion concentration is studied by recording the luminescence from various 
manifolds under the same experimental conditions. Comparison of the spectra allows for a 
determination of how an increase or decrease in a given dopant ion concentration affects 
the luminescence of the system as a whole, and helps to identify channels of excitation. 
The response of the system to changes in the power of the excitation source is also 
studied. The power of the excitation source affects the distribution of energy in the 
excited states of the system. Changing the power may increase or decrease the number of 


5 



excitations in a set of energy levels. Such changes represent a loss or a gain in the medium. 
Upconversion, for example, is a process in which an excited ion transfers its energy to 
another excited ion, exciting it to a higher level. This process results in a loss for one set of 
levels and a gain for another set of levels. Therefore, studying the response under changes 
in excitation power helps to identify which levels show a gain or a loss, and indicate how 
the levels are populated. 

The temporal response of the system under pulsed excitation at several different 
wavelengths has also been studied. The rate of change of populations in sets of energy 
levels has been studied with models that describe the observed behavior. Direct energy 
transfer from Tm -> Ho is accounted for in an analysis of rate equations. Cross-relaxation 
is accounted for in a model that incorporates spatial averaging over donor environments 
and diffusion among thulium ions. A descriptive analysis of the decay data is also 
employed in conjunction with other spectroscopic experiments to identify the modes and 
channels of excitation. 

1-5 QalLimt af Contents 

In the next chapter the theory of rare earth ions in solids is presented. Electronic 
configurations, the crystal field, energy levels, and the electron-phonon interaction is 
discussed. An outline of the principles involving the Judd-Ofelt theory are presented and 
methods of calculating emission cross sections are given. Finally, some models of energy 
transfer and their applicability are discussed. 

In chapter 3, YLF as a host lattice is presented. The structure and symmetiy of 
YLF crystals are discussed. The mechanical, thermal and optical properties have been 
collected from the literature and are tabulated here. 
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In chapter 4, the experimental apparatus and spectroscopic samples used in this 
study are presented. The details and specifications of absorption, cw emission and pulsed 
emission measurement apparatus are given in detail. Methods of measuring the polarized 
spectra are also given to unambiguously establish a proper reference for the spectra 
presented in this thesis. Finally the spectroscopic samples used in this thesis are listed, 
along with their RE ion dopant concentrations and growth method. 

In chapter 5, the results and interpretation for YLF:Ho are presented. The Judd- 
Ofelt theory is applied to YLFiHo and radiative rates are calculated for a number of 
levels. Absorption spectra and cw emission spectra in terms of cross sections are given. 
Finally, the temporal behavior of YLF.Ho to pulsed excitation is presented. 

In chapter 6, the results and interpretation for YLF:Tm are presented. The Judd- 
Ofelt theory is applied to YLF:Tm and radiative rates are calculated for a number of 
levels. Absorption spectra and cw emission spectra in terms of cross sections are given. 
Finally, the temporal behavior of YLF:Tm to pulsed excitation is presented, with special 
attention to characterizing and modeling cross-relaxation. Rates for the processes are 
calculated based on the model parameters. 

In chapter 7, the results and interpretation of co-doped YLF:Tm,Ho are 
presented. Sensitized luminescence is discussed and a descriptive analysis of the 
spectroscopic data is made. Upconversion is considered and rates are extracted from 
power dependent luminescence measurements. Decay data is presented and rate equation 
modeling is used to describe the Tm -> Ho energy transfer. 

The thesis ends with a summary of results and conclusions. Some suggestions for 
future work are also put forth. 
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CHAPTER 2 

Theory of Rare Earth Ions In Solids 


There is a theory which states that if ever anyone discovers exactly what the 
Universe is for and why it is here, it will instantly dissapear and be replaced 
by something even more bizarre and inexplicable. 

Douglas Adams 

The Restaurant at the end of the Universe 


2.1 Electronic Configuration 

The rare earths are divided into two sets of elements, the lanthanides and the 
actinides. The neutral forms of these elements have an electronic configuration made up 
by progressively filling the 4f electronic shell for the lanthanides and the 5f electronic 
shell for the actinides. A feature common to all the elements in the lanthanide series is a 
Xenon based electronic structure (ls 2 2s 2 2p 6 3s 2 3p 6 3d 10 4s 2 4p 6 4d 10 5s 2 5p 6 )with two or 
three outer electrons (6s 2 or 5d6s 2 ). The energy and spatial extent of the 4f wavefunctions 
show a sudden decrease at the beginning of the lanthanide series. Continuing through the 
lanthanide series, each additional added electron is accompanied by an increase in the 
effective nuclear charge. This causes a reduction of the 4f shell because of the imperfect 
shielding of one 4f electron with another. Each element of the lanthanide series is 
characterized by the number of 4f electrons it possesses, n=l for Cerium through n=14 
for Lutecium. The lanthanides are ionized by successive removal of electrons. The first 
stage of ionization consists of removal of a 6s electron, the second by removal of an 
additional 6s electron, and the third by removal of all the 6s and 5d electrons or perhaps a 
4f electron. This third stage of ionization comprise the triply ionized lanthanides, which 
are the most common valence states of RE ions in solids. 
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Many of the properties of rare earth ions result from the fact that the contracted 4f 
shell behaves as an inner shell. Since it behaves as an inner shell, the 4f wavefunctions do 
not penetrate into the environment outside the electronic configuration. The result is that 
the 5s 2 5 p 6 shell shields the f-electrons from interactions originating in the ions 
environment. 

2.2 IM£ryslal Field 

To study the energy level structure of electron configurations of rare earth ions in 
solids, the effect of the crystal lattice must be taken into account. When an ion is placed 
in a crystal lattice the free-space electronic energy levels are modified by interactions with 
the new environment. Typical absorption and emission spectra of RE ions in solids show 
groups of sharp lines indicating that these various groups of lines are due to transitions 
between 4f free ion states that are split by interactions with the surrounding charges of 
the lattice. In other words, the cumulative effects of the lattice charges surrounding the ion 
cause a perturbation that splits the ionic energy levels into sub levels. By splitting, it is 
meant that the 2J+1 fold degeneracy is partially removed. This is the well-known Stark 
effect. The perturbation that causes the splitting is an electrostatic field known as the 
crystal field. The theory that describes this effect is known as crystal field theory. In 
reality, however, a more general ligand field theory gives a more accurate picture. In this 
theory, covalent bonding plays an additional role since the electrons of the ion are 
partially shared with the surrounding ions or ligands. As was mentioned earlier, however, 
the shielding of 4f electrons from the surrounding environment reduces this effect. Ligand 
field theory is more appropriate to systems of transition metal ions, which are not 
shielded from the suirounding environment. Crystal field theory is then an accurate 
approximation when applied to RE ions. 
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If the above suppositions are correct then we must reconcile the above with the fact 
that transitions between 4f Stark levels are forbidden for electric dipole radiation by the 
parity selection rule, which states that electric dipole transitions can only occur between 
levels of opposite parity. Since the energy levels of the 4f shell all have the same parity 
there must be a mechanism that changes the parity of some of the levels. A resolution to 
this problem was found when it was recognized by VanVleck [Va37] that the crystal 
field, apart from splitting the energy levels by removing the 2J+1 degeneracy could also 
mix higher states of opposite parity into the 4f states. Odd power terms in the expansion 
of the crystal field potential can mix states of opposite parity into 4fi* from higher 
configurations, such as 4f N ' 1 5d, that make some of the transitions allowed. It should be 
noted that magnetic dipole transitions do not violate the parity selection rule (e g., they 
can occur between states of the same parity). A weaker mechanism that can also 
contribute to the admixing of states is the interaction of the electronic states with lattice 
vibrations. This mechanism becomes the dominant one when ions are placed in sites with 
inversion symmetry, since the crystal field contribution becomes zero in this case. 

2.3 Energy Levels of Rare Earth Ions in Solids 


The next question to be dealt with is to what extent the crystal field splits the levels 
and how many levels are obtained for a given crystal field. Wyboume [Wy65] gives an 
extensive account of this subject. 

In general, one starts with a free-ion Hamiltonian Hp given by : 


2 N N 2 N 2 

i=l i=l 1 i<j 1J 


(2-1) 


12 


where the first term is the sum of the kinetic energies of all the electrons of a 4f ion, the 
second term is the potential energy of all the electrons in the field of the nucleus, the third 
term is the repulsive Coulomb potential of the interactions between pairs of electrons, 
and the last term is the spin-orbit Hamiltonian, which accounts for coupling between the 
spin angular momentum and the orbital angular momentum. 

In the free atom there is a spherical symmetry and each level is 2J+1 degenerate. 
When the ion is placed in a crystal environment the spherical symmetry is destroyed and 
each level splits under the influence of the crystal field. In fact, the spherical symmetry is 
reduced to the point symmetiy at the ion site. The degree to which the 2J+1 degeneracy 
is removed will depend on the point symmetry surrounding the ion. This aspect will 
become clear shortly. The perturbed Hamiltonian for an ion in a crystal is written as: 

H = H f + V (2-2) 

where V, the perturbation Hamiltonian, is due to the potential provided by the crystal 
environment around the ion. Since the eigenfunctions of the free Hamiltonian possess 
complete spherical symmetry and are expressible in terms of spherical harmonics, it is 
natural to expand V in terms of spherical harmonics or operators that transform like 
spherical harmonics. So, the potential is expanded in terms of the tensor operators C kq : 

V = X B kq C kq(«i.<Pi) (2-3) 

k.q.i 

where the summation over i involves all electrons of the ion of interest. The terms B kq are 
the coefficients of expansion and are known as crystal field parameters. They are 
generally determined empirically from experimental data. It turns out that for identical 4f 
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electrons the energy matrix element of V between states of opposite parity (the crystal 
field split energies) is nonzero only for k even and k £ 6. This follows from the triangle 
rule for addition of angular momentum l j + 1 2 £ k. The values for k and q are also limited 
by the point symmetry of the ion. This arises from the fact that the Hamiltonian must be 
invariant under operations of the point symmetry group. Thus, the crystal field must also 
exhibit the same symmetry as the point symmetry of the ion, since it is part of the total 
Hamiltonian. Equating the crystal field expansion with the expansion that has been 
transformed through operations of the point symmetry group gives the crystal field for a 
particular ionic point symmetry. Thus, the spherical symmetry of an ion in a crystal is 
reduced to the point symmetry at the site of the ion. 

Once the point symmetry and the appropriate form of the crystal field are known, 
it is possible to construct the crystal field energy matrix. This matrix is then diagonalized 
using an estimated set of starting parameters. The resulting theoretical set of energy levels 
are compared to the set of experimental levels and, by an iterative fitting procedure, the 
parameters are adjusted to obtain the best overall fit to experiment. 

The thirty-two crystallographic point groups can be divided into four general 
symmetry classes as follows [Ru56]: 

1. Cubic: Oh , O , T,j,T h , T 

2. Hexagonal: D 6h , D 6 ,C 6v , Cgh , C 6 , D 3h , C 3h , D 3d , D 3 , C 3v , S 6 , C 3 

3. Tetragonal: D 4h , D 4 , C 4v , C 4 , D 2d , S 4 

4. Lower symmetry: D 2h , D 2 , C 2v , , C 2 , C s , S 2 , C t 

It can be shown that knowledge of the symmetry class at the site of the ion can be 
used to predict the number of levels a given J state splits into. The number of levels a 
state of angular momentum J splits into for a given point symmetry is presented in table 
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2-1. Alternatively, rare-earth ions can be used to probe the crystal symmetry if the 
number of levels of the ion can be determined. 


Table 2-1 Splitting Under Point Symmetry For Integral J 


£ 

Cubic 

Hexagonal 
Tetragonal 
Lower symmetry 


0 

1 

2 

3 


5 

6 

7 

_8 

1 

1 

2 

3 

4 

4 

6 

6 

7 

1 

2 

3 

5 

6 

7 

9 

10 

11 

1 

2 

4 

5 

7 

8 

10 

11 

13 

1 

3 

5 

7 

9 

11 

13 

15 

17 


As has already been stated, the positions of the levels arise from a combination of 
the Coulomb, spin-orbit, and crystal field interactions. These splittings of the 4f^ 
configuration are shown in figure 2-1. The electrostatic interaction leads to a split 2S+1 L 
on the order of 10 4 cm -1 . The spin-orbit interaction splits the levels further into 2S+1 Lj, 
separating the J states by 10 3 cm- 1 . Finally, the crystal field removes or partially removes 
the degeneracy in J yielding an energy level separation on the order of 10 2 cm* 1 . The 
extent to which the Stark split sublevels spread depends on the strength of the crystal 
field. The larger the crystal field, the larger will be the spread of the J sublevels. Figure 2-2 
is a more explicit example, Holmium LiYF 4 in this case, of how the states split under 
various interactions. 
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Figure 2-1 Energy Level Splittings Under Various Interactions 
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Figure 2-2 Energy Level Scheme ofHo 3 * LiYF 4 
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2.4 Electron-Phonon Inte r action 


When the lattice and the ion are treated as a coupled system, optical transitions are 
considered to occur between vibrational-electronic or vibronic states. Due to this extra 
dynamical interaction, the Hamiltonian in eq. 2-3 must be modified as: 

H = Hj + H P + Hjp (2-4) 

where Hj =H F + V is the energy of the ion (free ion plus static crystal field), H P is the 
phonon field energy, and H IP is the interaction energy between the ion and the lattice. 
This last term is the one responsible for transitions between crystal field levels through 
creation or annihilation of phonons. 

The crystal lattice can vibrate over a range of frequencies, producing relaxation by 
emission or absorption of phonons in that range of frequencies. These phonons usually 
have characteristic frequencies less than the energy difference between which the 
transitions occur, the J states. Several phonons are then required to conserve energy. This 
multiphonon relaxation is not the only manifestation of the oscillating crystal field. Line 
broadening, the appearance of vibrational sidebands in the optical spectra, the dependence 
of Hnewidth and the position of spectral lines on temperature are some other effects 
which result from the interaction of the lattice with the ion. 

2.5 Intensities of Rare Earth Spectra 

The optical spectra of rare earth ions in solids arise primarily from electric and 
magnetic dipole transitions between Stark split multiplets of the 4^ electronic 
configuration. Electric dipole transitions are parity forbidden between levels of the 4f** 
configuration, while those of magnetic dipole origin are allowed. It is known from 
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experiment, however, that the significant contributions to the intensities of most 
transitions are electric dipole in nature. This is due to the mixing of opposite parity, 
higher level energy configurations with the 4f** configuration, such as 4f N ' 1 5d. This 
admixing of configurations is caused primarily by the odd-order terms in the crystal Stark 
field. It should be noted, however, that even though magnetic dipole radiation is about six 
orders of magnitude smaller than electric dipole radiation in free ions, the electric dipole 
radiation occurs as a consequence of a perturbation for ions in a crystal. As such, some 
transitions show a significant magnetic dipole contribution. 

Judd [Ju62] and Ofelt [Of62] have shown that the line strength from the initial J 
manifold (\|fJl and the terminal J' manifold I xj/'J') can be expressed as a simple sum of the 
elements of tensor operators and phenomenological Judd-Qfelt parameters, which is 
given in eq. 2-5. 

Sed = InJ(vJlu ft V , ')| 2 (2-5) 

X= 2,4,6 1 1 


where <vJ I i U a) 1 1 vT> are the doubly reduced matrix elements. 

To arrive at this result, let 1 yJ) and I y'J') be two states of the 4fN configuration 
that become slightly admixed with states a = (y"J") from configurations like 4f N ‘ 1 5d. 
From perturbation theory we get two mixed parity states and *P b given by: 


('TaH'I'Jl- 


y (yJ| V|ct)(q| 
•^E( V J)-E(a) 


(2-6) 


l'5'b)=l¥'J')-S 

a 


|a)(a|V|yT) 

E(\|/T)-E(a) 


(2-7) 
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where V is the operator of the crystal interaction energy. The electric dipole operator can 
be written in a tensor operator form as: 

Pp = -e£ r i( c p), (2-8) 

i 

where p = 0 corresponds to the z-components for absorption or emission of ic-polarized 
light, and p - ± 1 gives the x ± iy components for o-polarized light. The matrix element of 
the electric dipole operator to first order in the crystal field perturbation is written as: 


(w Ip ^ <V J |Vla>(a|Pp|yT) ^ (yT|P p |aXa|V|y'J / ) 

\ al p| b/ Z, e(\j/J) - E(a) ^ E(y'J')-E(ci) 


(2-9) 


This equation does not lend itself to further calculations unless some 
approximations are made. Judd and Ofelt made two assumptions. The first was to assume 
that the states a=(\jr ,, j") are completely degenerate, which is equivalent to assigning to 
them an average energy E a = E(a). This allows one to remove the terms in the 
denominator from the sum and invoke the closure relation £ lot) (otl =1. If la) = |5d) 
then clearly this assumption constitutes a weak link in the theory since the energies of Sd 
states span a range of about 50,000 cm* 1 (from 50,000 cm* 1 to 100,000 cm -1 ) 
encompassing many different states with many different average energies. However, if the 
5d configurations are sufficiently complex and encompass many different terms, which 
when broken into individual sums over closely lying terms produced average energies that 
were roughly proportional to each other, then the error of this assumption would be 
reduced. Now, the second assumption is to take: 

E( ¥ J)-E 0 = E( V T)-E a (2-10) 
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This assumption is equivalent to the supposition that the 4f N_1 5d configurations lie 
far above the states involved in the transitions. This assumption is fairly well met, but 
carries with it the error associated with E a . At any rate, from these assumptions and the 
use of tensor operator algebra, one arrives at equation 2-5. 

The only other point to mention about eq. 2-5 is that the doubly reduced matrix 
elements vanish for X odd or X > 6. The values for the matrix elements can be found in 
various sources. The most complete set is given by Camall, Fields and Rajnak [Ca68a]. 
Experience has shown that the matrix elements vary only by small amounts from host to 
host, so they can be used regardless of the host materials studied. It should be noted that 
the expression for one-phonon vibronic transitions has an identical form as that of eq. 2-5 
and thus is incorporated into the Judd-Ofelt parameters. An equation similar to eq. 2-5 
can also be written for the magnetic dipole contribution to the line strength. 

s MD=(44HfKvI!L+2Slvl')l 2 (2-11) 


Equations 2-5 and 2-1 1 are related to the spontaneous emission probability. 


A(yJ; ¥ T) = 


64tt 4 e 2 
3h(2 J + 1)5? 



S E o + n S MD 


( 2 - 12 ) 


The connection between eq. 2-12 and experiment is provided by calculating the 
linestrength from absorption measurements. Specifically, the measured quantities are 
integrated absorption cross sections, that are related to the experimental linestrength by: 

J° a)<a (2 - b) 

manifold 
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The magnetic dipole contribution is given by the equation: 


S MD = 


3h(2J / + 1)X 
8jt 2 mc 


f MD 


(2-14) 


where is the magnetic dipole oscillator strength. These have already been calculated 
for many different transitions in rare-earth ions [Ca68b]. In cases where the magnetic 
dipole contributes significantly to the line strength it must be subtracted from the electric 
dipole part. 

To evaluate the spontaneous emission probabilities of eq. 2-12, the electric dipole 
and magnetic dipole line strengths must be calculated. This is easily done for transitions 
that terminate at the ground state from simple spectroscopy. The advantage of the theory 
of Judd and Ofelt is that once the set of intensity parameters {QjJ are found for a given 

ion-host combination, they can be used to calculate absorption and emission probabilities 
between any two fN levels. They are also useful for calculating radiative lifetimes for any 
excited manifold since this simply requires summing the transition probabilities of an 
excited state to every state below it down to the ground state. The way one determines 
the set of J-0 parameters is by performing a least squares fit between the theoretical 
expression in eq. 2-5 and the expression in eq. 2-13. More about this fitting procedure 
will be presented in a later chapter when it will be applied Tm 3+ and Ho 3+ in YLF. 


2.6 Emission Cross Sections 

The emission cross section at a given frequency can be determined from one of the 
forms of the Fiichtbauer-Ladenberg formula [Pa92]: 


o(v) 


X 2 

8jcn 2 X rad 


g(v) 


(2-15) 
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where g(v) is the spectral lineshape function, which can be written in terms of the 
intensity I(v) of an emission spectrum as: 


*(v) = j^ < 2 - 16 > 

It will be useful for purposes later to convert this into an equation involving the 
wavelength X. By noting that I(v) = I(X) and dv = (c / X )dX it can be shown that eq. 2- 
15 can be written as: 


o(X) = — — 


I(X) 


8jccn 2 T rad j^dX 


(2-17) 


It is also possible to obtain the emission cross section from a knowledge of the 
absorption lineshape and the energy levels of the ground state manifold and an excited 
state manifold. Equations connecting emission cross sections to absorption cross sections 
have been derived by McCumber[Mc64] and more recently by Payne et al. [Pa92], 

Emission and absorption cross section functions can be written as the sums of 
individual cross sections between the i th lower and j** 1 upper Stark levels multiplied by 
their respective probability of occupancy in thermal equilibrium at temperature T. 

^Em(v) = S f j°ji(v) (2- 18 a) 

ij 


°AB(v) = £ f i a ij( v ) (2- 18b) 

»j 


where, the Boltzmann fractions fi and fj can be written as: 
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(2- 19a) 


. _ gjg E|/kT _ gj -Ej/kT 
' ' Sg m e-EmAT - Z gnd 

m 


f 


j 


g j e~ Ej / kT 

Ig m e- Em/kT 


_ii_ e - E j/ kT 

^gnd 


(2-1 9b) 


where Zgnd, Z exc are the partition functions of the ground, excited states respectively. 

Figure 2-3 shows the Stark levels of a ground and excited state. A zero-line energy 
Ezl can be defined as the energy between the lowest Stark levels of the ground and excited 
manifolds. The energy difference between the lowest Stark level of the excited state and 
any given Stark level of the excited state can be written in terms of Ezl , Ei > and the 
transition energy hv. From figure 2-3, clearly: 

Ej=Ei+hv-EzL (2-20) 


Excited 

state 


Ground 

state 



Figure 2-3 Stark Levels of Ground and Excited state 
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Equation 2- 19b can now be written as 


f . = JLL e ~ E i/ kT e - (EzL ~ hv)/kT 

Z exc 


( 2 - 21 ) 


Comparing this with eq. 2-19a, eq. 2-21 can be written in terms of fi : 


f . - f. Ii.^El e -(EzL-hv)/kT 

1 j "M o. 7 c 
J ®i ^-exc 


( 2 - 22 ) 


Substituting this expression back into 2-1 8a, the emission cross section becomes: 

oem(v) = |i 0ji (v) (2-23) 

exc ij 1 


Invoking reciprocity ay =(gj/gi)Oji , the summation in eq. 2-23 becomes the 
absorption cross section. The final result is: 


OEM(v) = OAB(v)^e- (E ZL- hv >/ kT (2-24) 

^exc 


2.7 Energy Transfer 

When a solid-state material containing optically active ions is exposed to a radiation 
source, electronic states can be excited if the ions absorb some of the photons in the 
radiation field. At some later time the electronic states de-excite through the emission of 
photons, which show up as light, or they can give up phonons to the crystal lattice, 
which show up as heat. This simple picture works well if the ions can be considered as 
isolated systems, which would be the case for very small ion concentrations where the 
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ions are well separated from each other. As the concentration of dopant ions is increased, 
however, the spacing between ions becomes small enough to allow the ions to interact 
with each other. It is this interaction between ions that is responsible for the transfer of 
energy from one excited ion to another. Energy transfer processes can result from various 
mechanisms. Photoconductivity occurs when the absorption of photons in a material 
creates electron-hole pairs that can migrate in a solid carrying both energy and charge. 
This type of mechanism is found mostly in electrical conductivity. 

Another mechanism is radiative energy transfer in which a photon emitted by an 
excited ion is absorbed by another. This type of behavior is well understood in terms of 
the usual emission and absorption processes and will not be discussed here. However, the 
migration of radiative energy is important and does produce observable effects and cannot 
always be neglected. Generally this process can be taken into account by realizing that the 
radiative transfer of energy depends critically on the size of the solid-state sample and 
experimental geometry, so adjusting the experiment accordingly will minimize the effects. 

The type of energy transfer that this thesis will be concerned with is non-radiative 
energy transfer , in which the excitation energy from an ion of one species (the donor or 
sensitizer ) is transferred to an ion of another species (the acceptor or activator). This 
type of energy transfer can be viewed as a quantum mechanical resonance process 
involving the exchange of a virtual photon and can produce various effects on the 
luminescence of ions in solids (sensitized luminescence, fluorescence quenching, self- 
quenching, etc.) Non-radiative energy transfer can be a resonant one if the energy levels of 
the absorbing ions (donor) overlap the energy levels where emission occurs in the 
acceptors. The transfer can also be non-resonant if the phonons in the lattice are of 
sufficient energy to assist the transfer of energy by the creation or annihilation of 
phonons from the lattice such that the energy mismatch is compensated. 
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The mechanism of energy transfer between donor and acceptor was initially studied 
by Forster [Fo46,Fo48] and Dexter [De53] for the multipole and exchange interactions. 
The Hamiltonian describing the coupling between the electrons of the donor (D) and the 
electrons of the acceptor (A) due to the Coulomb interaction is: 



(2-25) 


where k is the dielectric constant and i and j are the electrons of the donor and acceptor 
respectively. This is the Van der Waals' interaction between ions in solution or solid. The 
initial and final states of the system can be described by a superposition of properly 
symmetrized single ion states. Forster calculated the transition rate using the Fermi's 
Golden Rule: 


WDA<E) = ^K¥ i |H DA |>l' t )| 2 p(E) <2-2«) 

The Coulomb Hamiltonian in eq. 2-25 can be expanded in a Taylor series about the 
vector representing the separation between the donor and acceptor. This leads to an 
expression for the transition rate as a sum of multipole interactions: 


pDD pDQ pQQ 

Wda =a“+ fi- + ii 77r + - 

A R 6 R 8 R 10 


(2-27) 


where C DD ,C°Q , and CQQ represents the strength of the dipole-dipole, dipole-quadrupole, 
and quadrupole-quadrupole interactions respectively. R is the distance separating the 
donor and acceptor ions. A large amount of experimental work based on the Forster- 
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Dexter model has been done in the past few decades. In particular, the Forster-Dexter 
theory has been used to determine which multipolar terms dominate in energy transfer 
and luminescence quenching processes, and to calculate values for the respective 
interaction parameters (C DD ,C°Q ,C*>, etc.). 

To extract relevant information about the microscopic interactions between acceptor 
and donor ions, a statistical analysis of the donor-acceptor distances in real, macroscopic 
crystals is necessary to account for the distribution of environments seen by different 
donor ions. Such an approach was originally taken by Forster [Fo49] and later generalized 
by Inokuti and Hirayama [InHi65], The general expression of the donor decay for any 
order of multipole interaction derived in [InHi65] can be written as: 


no(0 = n 0 (t)exp 


-:H* r H)N A ( c DA<) 


,3/q 


(2-28) 


where x is the lifetime of the donor in the absence of transfer, T(x) is the gamma function, 
Na is the acceptor ion density, Cda is an "interaction microparameter" for the donor to 
acceptor transfer, and q is the order of the multipole process (q=6 for dipole-dipole, q=8 
for dipole-quadrupole, and q=10 for quadrupole-quadrupole). In practice, this relation can 
be used to determine the order of the multipole process by fitting it with the experimental 
data. In reality, however, the similarity of the decay functions for q = 6, 8, 10 and the 
uncertainty inherent in the experimental data makes this approach inconclusive. Unless it 
is known that electric dipole transitions are forbidden, experimental decay data are usually 
fit with q=6 assuming a prevalent dipole-dipole process. 

There have been some objections to the applicability of eq. 2-28 in the 
interpretation of energy transfer. Grant [Gr71], in a paper that approached the 
macroscopic case of energy transfer with the use of rate equations, showed that the 
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concentration dependence of energy transfer is a result of the assumption that transfer 
occurs primarily between nearest neighbors. Grant points out that it is not clear why this 
is justifiable in principle, but probably does not introduce much of an error for small 
concentrations. Fong and Diestler [FoDi71], on statistical considerations, came to the 
conclusion that the dependence of transfer probability on concentration does not reflect 
multipole behavior, but is due to many body processes that reflect the number of ions 
participating in the process. More recently, Rotman [Ro89] has shown evidence for the 
correlated placement of ions in Nd:YAG, where energy transfer takes place between 
nonuniformly distributed ions in which a donor has an enhanced or excluded number of 
acceptors to interact with. 

Despite the objections made to it, eq. 2-28 has been found to describe the decay of 
donor ions quite well in many systems with low donor ion concentration. When the donor 
ion concentration is increased, the distance between donor ions decreases and it becomes 
more likely for excitation to migrate or hop to a near neighbor. Excitation energy may then 
be able to migrate among the donor ions before it is passed to an acceptor ion. With 
diffusion, excitation energy is transferred to an acceptor ion more efficiently by diffusing 
through the donors to interact with an acceptor at a much shorter distance, rather than 
interacting over a larger distance. The process of diffusion also leads to spatial equilibrium 
among the donor ions by averaging out the variation in environment as seen by the 
donors. For high concentrations of donors, the diffusion rate will be much faster than the 
D->A transfer rate and the average donor-acceptor interaction distance will be roughly the 
same over the entire system. This results in D->A transfer rates that are effectively the 
same over the entire ensemble of ions. The rate limiting step for the donor relaxation then 
becomes either the donor-acceptor transfer rate or the acceptor relaxation rate. The 
outcome is an exponential decay rather than the t 1/2 exhibited by eq. 2-28. This is the case 
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for fast diffusion and is ideal for a treatment with rate equations where the spatial 
variations in donor environments are not an inherent property of rate equations. 

Experimentally, it is very common to see decays that exhibit non-exponential 
behavior at early times and then a longer exponential tail. For moderate concentrations the 
rate of energy diffusion within the donor system may not be sufficient to maintain a 
uniform excitation of donors. In this case, certain donor sites will interact with acceptors 
over a larger distance than others that see enhanced rates of interaction due to diffusion. 
The decay of the donor system is then composed of competing processes. This is the 
case of diffusion-limited relaxation and has been treated by Yokota and Tanimoto 
[YoTa65]. They solved the diffusion equation of the form: 


f=- Dv: 


C DA 


n _£_ n £±M 

x j |r-rj 


(2-29) 


where the first term is the rate of change due to diffusion with diffusion constant D, the 
second term is the natural lifetime (radiative + nonradiative) in the absence of cross 
relaxation and diffusion, and the third term is the interaction with i acceptors at position 
ri, with the dipole-dipole interaction parameter Cda- Using the method of Pade 
approximants they found the following solution: 


n(t) = n(0)exp 



ht H2( l+10.87x+15.5x 2 f /4 
Dt ( l+8.743x ) 


(2-30) 


where b = ^ 7t 3 ^ 2 N ^ ^/C DA t and x = DC D ^ 3 t 2 ^ 3 


When D->0, eq. 2-30 reduces to eq. 2-28 for a dipole-dipole case (q = 6). 
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At early times, t«D _3/2 CDA 1/2 » the diffusion is unimportant and the decay is non- 
exponential as in eq. 2-28. As t-><», the donors will decay exponentially at a rate 
determined by the migration. As the migration becomes faster the boundary between 
these two types of behavior shifts to shorter times until the decay appears to be purely 
exponential. This changing time dependent decay is a distinguishing characteristic of 
diffusion limited relaxation. At veiy long times, the fluorescence decay can be written as 
[YoTa67]: 

n(t) = n(0)exp[-±--J-j (2-31) 

where, 

^ = 0.5l(4jtN A C^D 3/4 ) (2-32) 

So, at later times in the decay stage, the donors with accessible acceptor 
environments have all transferred their energy to the acceptors. Further decay is limited 
until migration brings the donor excitation close enough to some acceptors to facilitate 
transfer. The decay is then seen to be limited by the diffusion rate. It is exponential and 
proceeds at the rate given by eq. 2-32. This regime of the decay is most easily 
investigated because analyzing decay curves at longer times is easier since they become 
exponential and generally comprises a substantial portion of the decay profile. 
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CHAPTER 3 

LiYF 4 (YLF) as Host Lattice 


God has put a secret art into the forces of Nature so as to enable 
it to fashion itself out of chaos into a perfect world system. 

Immanuel Kant 
Universal Natural History 
and Theory of the Heavens 


3.1 Stmtfure. flf Sffllidsi th e C rystal latt i ce 

A crystal can be described as a collection of atoms or molecules arranged in a 
periodic pattern. This periodic pattern extends over an interval that is large compared to 
the unit of periodicity, the so-called unit cell. The structure of any crystal can be 
described by defining a set of basis vectors along the crystallographic axes. In this way a 
lattice can be generated from a given point by repeated application of the three basis 
vectors, x with magnitude a, y with magnitude b, and z with magnitude c. The points of 
the lattice defined in this way do not have to correspond to actual positions of atoms, but 
are meant to indicate that the environment about any one lattice point is identical to that 
around any other lattice point. The position of the point lattice can be chosen so that each 
point corresponds with the positions of one of the kinds of atoms which make up the 
crystal. Once the placement and number of atoms within a unit cell are known, the entire 
structure of the crystal is determined by translational periodicity. 

The entire problem of determining the structure of a crystal is reduced to finding 
the positions of a relatively few number of atoms within a unit cell. The entire crystal is 
then visualized as made up of multiple unit cells, sharing edges, comers or faces. 
Crystallographic axes are chosen in a way that shows the symmetry of the lattice. 
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3-2 Symmetry in Crystal Systems 


By considering the underlying symmetry of crystal systems, the classification and 
detailed description of various types of crystals is greatly facilitated. The point lattice 
system described above serves as a basis of classification of crystals into one of seven 
systems (Triclinic, Monoclinic, Orthorhombic, Tetragonal, Cubic, Hexagonal, and 
Rhombohedral). This classification can be extended by considering unit cells that contain 
lattice points in addition to those at the corners. For example, consider a simple cube. 
This is the primitive cubic. A cube with a lattice point at its center is a body-centered 
cubic system, and a cube with a lattice point in the center of each face of the cube is a face 
centered cubic system. Each of these classifications is known as a Bravais lattice , of 
which there are fourteen. The reader is referred to a book on solid state physics 
[AsMe76] for more details. 

The point here is that an examination of the characteristics of well-formed crystals 
may reveal certain elements of symmetry that can be exploited with the powerful 
methods of group theory. An object is said to possess symmetry if an operation such as 
rotation, or reflection about a particular axis transforms the object into a configuration 
indistinguishable from the original. It turns out that there are 32 crystal classes that can be 
obtained by considering the symmetry operations that transform the crystal lattice to its 
original configuration with one of the lattice points remaining fixed. These are known as 
the 32 point groups since they refer to the symmetry about a single lattice site. Further, 
translatory movements of the lattice in a given direction and subsequent symmetry 
considerations of the new arrays obtained lead to a set of 230 space groups. All crystals 
must belong to one of the 230 space groups. Knowledge of the presence of symmetry 
within the unit cell, and X-ray diffraction methods are useful in determining the number of 
atoms and their relative positions within a given crystal lattice. 
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3.3 The Structure of UYFt 


Lithium Yttrium Fluoride, LiYF 4 (YLF) is a simple fluoride crystal with a 
tetragonal scheelite structure (space group C* A — / a, number 88) [Th61], [Sh69].The 

crystal structure of YLF is pictured in figure 3-1. 

Rare earth impurities introduced into YLF replace the Yttrium ions 
substitutionally. The point group symmetry at the site of the Yttrium ions is S 4 This 
means that while the overall symmetry of YLF is tetragonal, the symmetry at the site of 
the ions in the crystal is S 4 , so called because there exists a 4-fold rotation-reflection axis. 
A 4-fold rotation-reflection axis is one in which a rotation about it equal to n/2 followed 
by a reflection in the plane perpendicular to the axis produces no change. This seemingly 
technical detail is very important in determining selection rules for dipole-dipole 
transitions in laser crystals. 

It can be seen from figure 3-1 that the crystallographic directions in YLF are not 
equivalent. There exists two equivalent directions (the a-axes), and one independent of the 
other two (the c-axis). Crystals that exhibit this property are called anisotropic. YLF is an 
example of a uniaxial anisotropic crystal, as opposed to a biaxial crystal that has three 
independent directions. 

Trivalent rare-earth ions replace, substitutionally. Yttrium (Y 3+ ) ions in YLF with 
little change in lattice parameters since rare-earth ions and Yttrium have nearly the same 
atomic radii. The lattice constants for YLF have been measured from x-ray diffraction data 
[BlSa79], where they are reported to be a = 5. 166 A and c = 10.733 A. 

Using these lattice constants, table 3-1 shows the calculated distances from 
reference ion 0 (see figure 3-1) to neighboring Yttrium ions, and the corresponding number 
of nearest neighbors of types 1,2,3 and 4. 
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Table 3-1 Nearest Neighbor Distances in UYF 4 (YLF) 


HESEfflH 

Distance from Ion 0 (A) 


1 

3.7247 

i 4 

2 

5.1668 

4 

3 

6.3694 

8 

4 

6.9420 

8 

2 

7.3069 

4 

1 

8.2015 

4 

2 

8.4541 

4 

3 

9.6934 

8 

4 

9.7744 

16 

2 

9.9080 

8 
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Figure 3-1 Crystal Structure of UYF 4 (YLF) 



3.4 imtei af Refraction of LiYFi 


The index of refraction of YLF crystals has been measured previously [CaLi75] 
for both ordinary and extraordinary waves in the wavelength range 0.225pm to 2.6pm. 
The data reported in this reference is listed in table 3-2 and presented graphically in figure 
3-2. The data in table 3-2 has been fit by the method of least squares to the Sellmeier 
equation: 


n 2 = A + - 


BX 2 DX 2 


(X 2 -C) (Xf-E) 


(3-1) 


The Sellmeier coefficients were found to be: 
ordinary (n 0 ): A=1.35272, B=0.70757, C=0.00893, D=0.49161, E=124.992 

extraordinary (n e ). A=1 .35743, B=0.80225, C=0.00916, D=0.59253, E=145.892 


The units on C, and E are pm 2 , while A, B and D are unitless. The root mean 
square deviation for the fit of the refractive index data to eq. 3-1 was a - 0.000040. The 
temperature dependence of the refractive index of YLF in the visible region of the 
spectrum has also been done previously [BaGe80], where it was found that n decreases 
by small amounts at lower temperatures. The small negative variation of the refractive 
index tends to minimize thermal-focusing effects in solid-state laser rods, which can cause 
degradation of laser beam quality by reducing the TEMoo mode radius, allowing higher- 
order modes to oscillate. 

It should be pointed out that the terms ordinary and extraordinary index of 
refraction refer to the direction of the propagating electric field with respect to the c-axis 
of the crystal. Ordinary refers to an electric field perpendicular to the c-axis, and 
extraordinary refers to an electric field parallel to the c-axis. 
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Table 3-2 Refractive Index Values in LiYF 4 (YLF) 


Wavelength (|xm) 

n e 

n o 

2.600 

1.4602 

1.4381 

2.400 

1.4617 

1.4397 

2.200 

1.4632 

1.4411 

2.000 

1.4645 

1.4424 

1.800 

1.4658 

1.4437 

1.600 

1.4670 

1.4449 

1.400 

1.4681 

1.4460 

1.200 

1.4693 

1.4471 

1.000 

1.4708 

1.4485 

0.900 

1.4716 

1.4492 

0.800 

1.4726 

1.4502 

0.700 

1.4741 

1.4516 

0.600 

1.4962 

1.4502 

0.500 

1.4796 

1.4567 

0.450 

1.4822 

1.4590 

0.400 

1.4822 

1.4624 

0.350 

1.4913 

1.4674 

0.300 

1.5001 

1.4754 

0.250 

1.5156 

1.4895 

0.225 

1.5287 

1.5014 
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3.5 Mechanical Properties ofLiYF 1 


The mechanical properties of YLF have been taken from various sources in the 
literature. They are listed below in table 3-3 


Table 3-3 Mechanical Properties of UYF4 (YLF) 


Young's modulus 

B = 8.0xl0 10 N/m 2 

ref. [BlSa791 

Volume compressibility 

Pv = 12.5xl0* 12 N/m 2 

ref. [BlSa791 

Linear compressibility 

Pz= 3.4xl0* 12 N/m 2 
Pxy = 3.4xl0* 12 N/m 2 

ref. [BlSa79] 

Density 

p = 3.99 gm/cm 3 

ref. fBlSa791 

Poisson's ratio 

V = 0.33 

ref. [Ko92] 

Strength 

3.3x10*7 N/m 2 

ref. [Ko92] 

Hardness (ICnoop) 

260-325 kg/mm 2 

ref. [ChFo731 

Hardness (Mohs) 

4-5 

ref. [Ko92] 


3.6 Thermal Properties of LiYFj 

The thermal properties have also been taken from various sources in the literature. 
They appear in table 3-4. It should be mentioned that the phonon energies for YLF in 
table 3-4 represent the high frequency energies of the phonons of the effective modes 
involved in multiphonon nonradiative transitions. The range of phonons observed from 
Raman spectra [MiRa70], [CaRa75] is 140 - 560 cm* 1 , but due to their small mass the 
Lithium atoms are expected to vibrate at their highest frequencies for a given mode. The 
mode assignment labels for the phonons given in the above references are for C411 
symmetry not S4. The phonon modes and their corresponding highest frequency phonons 
are given in table 3-5. 
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Table 3-4 Thermal Properties of L 1 YF 4 (YLF) 


Thermal Conductivity 

k- 1.5xlO' 2 cal./cm-sec-°C 

BWHIBBBi 

Coeff. of Linear expansion 
(0-100°C), AT=100°C 

a-axis: a= 13.8x10*6 °C _1 
c-axis: a = 9.0x10*6 ’C * 1 

ref. [ChFo73] 

Melting Point 

T m = 825 °C 

1 

Debye Energy 

0D = 279cm* 1 



Ep = 424 - 566 cm * 1 

ref. [MiRa70] 

dno/dt 

-0.54x1 O ' 6 (@435.8 nm) 
-0.67x1 0* 6 (@546.1 nm) 
-0.91x10*6 (@578.0 nm) 

ref. [BaGe80] 

dne/dt 

-2.44x10*6 (@435.8 nm) 
-2.30x10-6 (@546. lnm) 
-2.86x10-6 (@578.0 nm) 

ref. [BaGe80] 


Table 3-5 Phonon Modes and Phonon Energies in UYF 4 (YLF) 


Mode 

Highest frequency phonon 

A e 

425 cm * 1 

Bg 

427 cm - 1 

Eg 

446 cm * 1 

A u 

490 cm - 1 

Eu 

566 cm * 1 
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Figure 3-3 Linear Thermal Expansion of YLF 


Figure 3-3 shows a graph of the linear thermal expansion of YLF. The coefficient 
of linear expansion can be calculated from the data in table 3-3. It is defined as: 


1 AL 
ATL 0 


(3-2) 


The coefficient of linear expansion for the a-axis and c-axis in YLF appears in table 
3-4 for a temperature range AT = 100 °C. 
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CHAPTER 4 

Experimental Apparatus and Samples 


In considering the intent for which a clock has been designed, I can scry 

that it is false to its nature when it does not correctly indicate the time. 

Rene Descartes 

Meditations on First Philosophy 

4.1 Absorption Measurements 

Absorption measurements were made with a Perkin-Elmer Lambda9 UV/Vis/NIR 
Spectrophotometer. The Lambda9 incorporates two juxtaposed monochromators. Each 
one of them is equipped with two large gratings, a photomultiplier and, a PbS (lead 
sulfide) detector. The range of wavelengths extended from 190-3200 nm. The schematic of 
the optical system is shown in figure 4-1. 

The excitation sources included a deuterium lamp for the UV range and a 
Tungsten-Halogen lamp for the visual and infrared regions. The instrument resolution was 
.05 nm due to limitations on the data acquisition software. Careful selection of various 
control parameters on the Lambda9 allowed a compromise to be made between spectral 
resolution and signal to noise ratio. 

A double depolarizer/polarizer assembly was used since the crystals analyzed in 
this study were uniaxial. Uniaxial crystals exhibit varying degrees of transmission for 
different orientations and polarizations. The double depolarizer/polarizer assembly is a 
complete sample compartment that fits into the Lambda9 as one unit. The depolarizer 
mounts in the optical path between the entrance slit and the crystal. This depolarizer was 
necessary since light dispersed from a grating is partially polarized. The depolarization 
effect, however, is not a true depolarization, but a statistical mixing of all polarizations 
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Figure 4-1 Schematic of Absorption Spectrophotometer Optical System 
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within the cross-section of the radiation beam. The double polarizer mounts between the 
crystal and the exit slit leading to the detector. It consists of two wedge-shaped Glan- 
Thompson prisms that can be rotated between 0 and 180 degrees by a mechanical drive 
located outside the sample compartment. 

The crystals were placed in the sample compartment with the c-axis of the crystal 
oriented in a vertical direction perpendicular to the optical path of the excitation beam. By 
orienting the mechanical drive to 0°, the tc-polarized (electric field of the optical beam 
parallel to the c-axis) transmission spectra was measured. By orienting the mechanical 
drive to 90° the o-polarized (electric field of the optical beam perpendicular to the c-axis) 
transmission spectra was measured. See figure 4-2 for orientation configurations. 

Due to the optical materials used in the depolarizer/polarizer the measurable 
wavelength range was 300-2600 nm. The transmission spectra were corrected for Fresnel 
losses at the crystal faces and converted to absorption spectra. 

4.2 Continuous Luminescence Measurements 

The apparatus used for continuous luminescence measurements consisted of an 
excitation system, a collection system, a detection/ amplification system, and a 
processing/recording system (see figure 4-3 for details). 

4.2.1 The Excitation System 

The excitation source used was a SONY model SLD 303XT-1 GaAlAs double- 
hetero laser diode. The laser diode assembly consisted of the photodiode, which had a 
thin active GaAlAs layer .05pm thick, a built in TE cooler for wavelength tuning, and a 
thermistor for temperature monitoring. The current for the lasers operation was provided 
by a Spectra Diode Labs SDL800 Laser Driver. The SDL800 had a current controller that 
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In the absorption measurements , 
the optical beam is depolarized 
before passing through the cry stal 
resulting in a randomly directed 
E-field vector. Polarization select 
is made by orienting a polarizer ii 
the desired direction after the bea 
is attenuated upon passage throug 
the crystal. See below. 




After the beam passes through 
a horizontal polarizer, the 
vertical component of the E-field 
is eliminated. The measured 
field is denoted as a-polarized. 
(Note that E is perpendicular to 
the c-axis.) 




After the beam passes through a 
vertical polarizer, the horizontal 
component of the E-field is 
eliminated. The measured field i: 
denoted as 7 t-polarized. 

(Note that E is parallel to the 
c-axis.) 


Figure 4-2 Polarized Absorption Measurements 
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Figure 4-3 Continuous Luminescence Apparatus 
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was limited to 1 ampere, and a temperature controller to tune the wavelength in the range 
785 ±15 nm. The current limit was set to 980 mA and the setpoint temperature to 6.8°C 
for an operating wavelength of 785.5 nm and an output power of 350m W. 

Laser diode beams are highly divergent, making it difficult to collect and collimate 
the beam. To overcome this an 8mm Melles-Griot lens was used in conjunction with an 
anamorphic prism pair to collect and collimate the laser radiation into a circular region 
with a 5mm diameter. The loss of diode light was kept at a minimum by using coated 
optics and by careful positioning of the lens/prism unit. This circular beam was then 
reflected from a mirror and focused onto the crystal with a 5 cm plano-convex calcium 
fluoride lens. This produced an excitation along the c-axis of the crystal after focus to a 
spot size measured to be 140pm x 260pm in diameter. 

4.2.2 Optical Collection System 

The luminescence from the crystal was collected by a 6 cm diameter (f=7cm) 
convex lens placed at a distance of 7 cm from the crystal to form a collimated beam, which 
was then focused by a 25 cm focal length lens onto the entrance slit of the spectrometer. 

The spectrometer used for wavelength selection of the luminescence signal was a 
0.5 meter McPherson model 216.5. The gratings used were a 1200 grooves/mm blazed 
@500 nm for measurements in the visible (300-850 nm), and a 600 grooves/mm blazed 
@1.6pm for measurements in the IR (850-2100 nm). Power dependence of these 
measurements was done by placing a variable attenuator between the laser and reflecting 
mirror. Polarization measurements were done by placing a polarizer between the focusing 
lens and the entrance slit of the spectrometer. The entrance and exit slits were kept at 
50pm in the visible and 500pm in the infrared. Figure 4-4 illustrates the polarized 
emission measurements. 
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In the case of polarized emission studies pictured above, the laser pumps along the c-axis of 
the crystal .exciting optically active ions in the crystal. While the laser is a coherent source 
with a definite polarization, the luminescence from the crystal is not. The incoherent 
radiation produced by spontaneous emission is unpolarized and spreads in all directions. 
Polarization selection is made by orienting a polarizer in the path of the collected 
luminescencse. The top diagram is an example of Jt-pol arization (E is parallel to the c-axis of 
the crystal.) The bottom diagram is an example of a-polarization (E is perpendicular to the 
c-axis of the crystal.) It should be noted that the the important aspect is the direction of the 
polarized E-field, and the propagation vector k relative to the c-axis. 


Figure 4-4 Polarized Luminescence Measurements 
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4.2.3 Detection and Amplification System 


The detectors used at the exit slit of the spectrometer were a Perkin-Elmer 
photomultiplier for visible wavelengths and an Electro-Optical Systems liquid nitrogen 
cooled Indium Antimonide detector for IR wavelengths. 

Signal amplification was performed by a PAR (Princeton Applied Research) 
model HR-8 lock-in amplifier in conjunction with a PAR model 125 mechanical chopper. 
The high voltage supply used for the photomultiplier was an EMI Gencom model 3000R. 
A voltage of 1250V was used in this experiment. The IR detector was powered with two 
6V batteries. 

4.2.4 The Processing and Recording System 

The analog signal from the lock-in amplifier was processed through connection to 
a DSP Technologies Crate System. The signals from the lock-in and chopper were 
transformed by passing them through a CAMAC module 1030 Amplifier. This provided 
the necessary voltage couplings (AC or DC), voltage offsets and signal attenuation for 
passage to a CAMAC module 2012 transient digitizer. The 2012 module sampled the 
analog signal and converted it to a digital number by sampling a given number of pre- 
trigger signals. The digital signal was then sent from the CAMAC module 6001 crate 
controller through an RS-232 connection to a Gateway 2000 computer system. A 
computer program was written to display the output data. Controls of the scanning 
parameters were entered at the keyboard and sent to the crate controller for initialization. 

4.3 Pulsed Luminescence Measurements 

Two laser sources were used for pulsed excitation measurements. Both of these 
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lasers were pumped with an Nd:YAG laser. The first was a 532 nm Nd:YAG pumped 
Ti:Al203 (Titanium Sapphire). The optical schematic of this laser is shown in figure 4-5. 
The Ti:Al203 laser was tunable over the wavelength range of 600 - 900 nm. The tuning 
was achieved with a prism in the laser cavity. When the laser was tuned to 780 nm the 
observed pulse duration was <100 ns and the peak output power was 25 mJ. The second 
laser source was a 1.33 pm NdiYAG pumped CoMgF2 (Cobalt Magnesium Fluoride). 
The optical schematic of this laser is shown in figure 4-6. The CoMgF2 laser was tunable 
over the wavelength range 1.7 - 2.1 pm. The tuning was achieved with a Brewster's angle 
two-plate birefringent filter in the laser cavity. When the laser was tuned to 1 .75 pm the 
observed pulse duration was 240 ns at HWHM and the peak output power was 
approximately 6 mJ. The size of the beam was measured with a knife edge technique 
[Bi85]. The beam was found to have a spot size of 200 x 250 pm. 

Wavelength selection of the fluorescence was achieved with a 0.75m Spex 1702/04 
spectrometer equipped with 1200 grooves/mm grating for the visible region of the 
spectrum. Narrow band filters were used for wavelength selection in the IR spectral 
region. 

The luminescence signals after wavelength selection were monitored with a 
Hamamatsu GaAs photocathode photomultiplier in the visible region of the spectrum and 
an EG&G Optoelectronics Indium Antimonide detector in the IR region of the spectrum. 
The Hamamatsu PMT was housed in a Products for Research, Inc. housing that could be 
cooled with water to 40°C below ambient temperatures. The EG&G detector was liquid 
nitrogen cooled to 77 K°. 

The voltage signals from the detectors were recorded with a LeCroy CAMAC 
6810 waveform recorder. The Catalyst software package for the LeCroy CAMAC was 
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Figure 4-5 Pumping System and Optical Schematic of the Ti :Al 203 Laser 



Figure 4-6 Pumping System and Optical Schematic of the CoMgF 2 Laser 
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used to perform an on screen average of the digitized signal waveform. In this way, 
fluctuations from background noise could be eliminated by averaging over say 500 or 1000 
waveforms. No preamplifiers were used to boost the signal level and the inputs from the 
detectors were terminated with a 50Q resistor to eliminate distortions of the signal 
waveform by transients in the equipment circuitry. 

4.4 Temperature Controls 

The temperature controlling system used consisted of four interconnected units: 
an RMC Cryosystems coldhead, an RMC Cryosystems model EC-2 compressor, a 
Veeco VS-9 vacuum pumping station, and a Lakeshore 320 Autotuning Temperature 
Controller. The apparatus is shown in figure 4.5. 

The vacuum pumping station serves to remove air and other impurities from the 
coldfmger cavity to better than 10~ 3 Torr. This is necessary to minimize convective heat 
conduction and achieve cryogenic temperatures less than 70K°. 

The cryogenic system operates as a closed-cycle refrigeration unit using 
compressed helium gas that expands inside the coldhead. The heat is removed by passing 
the helium through a heat exchanger to cool the crystal below ambient temperatures. 

The temperature inside the coldfinger was monitored with a Lakeshore 320 
temperature controller, which contained a silicon diode. The Lakeshore unit was 
autotuning down to about 23°K where the accuracy was rated at 23°K ± ,5°K. The unit 
was able to control down to 10°K with a reduced accuracy of ± 1°K. 

4.5 SpgctroscQpi f S amples 

The spectroscopic samples used in this study were crystals of Lithium Yttrium 
Fluoride, LiYF 4 (YLF), doped with triply ionized rare earth ions of Thulium (Tm 3+ ) and 
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Holmium (Ho 3+ ). Both singly doped and co-doped samples were used. A list of the 
samples available for study and the corresponding percent concentrations of the dopant 
ions is given below. 


a. ) Tm(4%)Ho(0.5%) YLF 

b. ) Tm(4%)Ho(l .0%) YLF 

c. ) Tm(4%)Ho(l .5%) YLF 

d. ) Tm(5%)Ho(0.2%) YLF 

e. ) Tm(6%)Ho(l%) YLF 


f. ) Tm(12%)Ho(l%) YLF 

g. ) Tm(0.2%)Ho( 12%) YLF 

h. ) Tm(0.5%) YLF 

i. ) Tm(3%) YLF 

j. )Tm(12%) YLF 


Samples a, b, c, d, f, g, i } and j were acquired from Oklahoma State University and 
were grown by the Bridgman method. Sample e came from Optical Lightning Corporation 
(growth method unknown). Sample h was from M.I.T and was grown using the 
Czochralski method. All the samples were compared to each other for accuracy of their 
dopant concentrations. It was found that samples f, g, and j, were in error from what the 
manufacturer specified. They were compared to samples whose dopant concentrations 
were known accurately, and corrected accordingly. The corrected values are the same as 
those given above. Throughout this thesis samples a-k will be referred to without the 
percentage sign. For example T m(4%)Ho(0 .5%) will be referred to simply as 
Tm(4)Ho(.5). 


57 



Rfiferencss to Chapter . A 


[Bi85] H R Bilger and T. Habib, "Knife-edge scanning of an astigmatic Gaussian Beam," 
Appl. Optics 24, 686-690 (1985) 


58 



CHAPTER 5 

YLF: Ho - Results and Interpretation 


The views of space and time which I lay before you have sprung from 
the soil of experimental physics, and therin lies their strength. 

Hermann Minkowski 
Space and Time 


s.i Eii&rgy Levels of YLF; Ho 


The energy levels for Ho 3+ in LiYF 4 have previously been measured and are 
reported in [KaWo76], In this thesis the experimental set of energy levels found in this 
reference will be used exclusively. On the basis of these energy levels, an energy level 
diagram was constructed. This is pictured in figure 5-1. 
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Figure 5-1 Energy Level Diagram of YLF:Ho 
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Tables for the upper and lower transition wavelengths (table 5-la) and energies 
(table 5-lb) for a given manifold in YLF:Ho have been made from the energy levels given 
in [ KaWo76], These tables are useful when looking for a resonance between a given pair 
of manifolds. 


Table 5-la Range of Transitions (nm) in YLF: Ho 


£X£ 512 515 515 514 5F5 5S2 5F4 5F3 5F2 


518 

u 

0 



1889 

WAV 

1136 

VAW 

882 

VLA-3 

738 

UL W 

638 

WWA 

539 

WiL_3 

534 

481 

471 


0 

2066 

1196 

915 

777 

659 

550 

546 

492 

480 

517 

- 

0 

2744 

1618 

1193 

952 

747 

737 

640 

622 


- 

0 

2961 

1681 

1267 

981 

758 

751 

652 

631 

516 

- 

- 

0 

3755 

2056 

1431 

1014 

996 

827 

797 


- 

- 

0 

4091 

2278 

1494 

1032 

1020 

845 

811 

515 

- 

“ 

- 

0 

4363 

2265 

1372 

1339 

1050 

1002 


- 

- 

- 

0 

5402 

2408 

1399 

1376 

1075 

1021 

514 

- 

- 

■- 

- 

G 

4048 

1872 

1810 

1320 

1245 


- 

- 

- 

- 

0 

5120 

2021 

1973 

1409 

1318 

5F5 

- 

- 

- 


- 

0 

3290 

3103 

1897 

1745 


- 

- 

- 

“ 

- 

0 

3541 

3395 

2010 

1830 

5S2 

- 

- 

- 

- 

- 

- 

0 

43668 

4391 

3656 


- 

- 

- 

- 

- 

- 

0 

133333 

4750 

3855 

5F4 

- 

- 

- 

- 

- 

- 

- 

0 

4638 

3825 


- 

- 

- 

- 

- 

- 

- 

0 

5202 

4147 

5F3 

- 

- 

- 

- 

- 

- 

- 

- 

0 

17123 


- 

- 

- 

- 

- 

- 

- 

- 

0 

27548 



-518 

. 517 

Table 5- lb. Range of Transitions (car 
5X6 515 5X4 5F5 

») of YLF: Ho 

_ 5F3 

5F2. . 

518 

0 

5293 

8796 

11333 

13532 

15654 

18524 

18707 

20755 

21213 


0 

4838 

8355 

10925 

12869 

15170 

18163 

18284 

20314 

20803 

517 

- 

0 

3643 

6180 

8379 

10501 

13371 

13554 

15602 

16060 


- 

0 

3377 

5947 

7891 

10192 

13185 

13306 

15336 

15825 

516 

- 

- 

0 

2663 

4862 

6984 

9854 

10037 

12085 

12543 


- 

- 

0 

2444 

4388 

6689 

9682 

9803 

11833 

12322 

515 

- 

- 

- 

0 

2292 

4414 

7284 

7467 

9515 

9973 


- 

- 

- 

0 

1851 

4152 

7145 

7266 

9296 

9785 

514 

- 

- 

- 

- 

0 

2470 

5340 

5523 

7571 

8029 


- 

- 

- 

- 

0 

1953 

4946 

5067 

7097 

7586 

5F5 

- 

- 

- 

- 

- 

0 

3039 

3222 

5270 

5728 


- 

- 

- 

- 


0 

2824 

2945 

4975 

5464 

5S2 

- 

- 

- 

- 

- 

- 

0 

229 

2277 

2735 


- 

- 

- 

- 

- 

- 

0 

75 

2105 

2594 

5F4 

“ 

- 

- 

- 

- 

- 

- 

0 

2156 

2614 


- 

- 

- 

- 

- 

- 

- 

0 

1922 

2411 

5F3 ' 

-- 

- 

- 

- 

- 

- 

- 

- 

0 

584 


- 

- 

- 

- 

- 

- 

- 

- 

0 
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5-2 JyddzQfelt Analysis o fY LFiffo 


The Judd-Ofelt theory has become a standard tool for the calculation of radiative 
transition probabilities of rare earth ions in solids. The basic principals and relevant 
equations of this theory have already been given in chapter 2. The application of this 
theory to YLF.Ho will be presented in this section. 

Absorption measurements were made on a sample of YLF doped with 0.5% 
Holmium - Ho(.5). Since YLF is an anisotropic crystal that exhibits different absorption 
intensities for different polarizations, two sets of spectra were collected over the 
wavelength range 300 - 2300 nm. The absorption spectra for Jt-polarization appear in 
figures 5-2, and the absorption spectra for o-polarization appear in figures 5-3. 

The measured quantities in the Judd-Ofelt theory are the integrated absorption 
profiles over a given manifold, but since there are two polarizations to consider, they 
must be suitably averaged to find the total integrated absorption. Toward this end, a 
modified Beer-Lambert law for absorption attenuation can be written as (Appendix A): 

I(X) = I 0 {^exp[-a Jc a)N^+|expC-a 0 (X)N£]} (5-1) 

A spatially averaged Beer-Lambert law can be written as I(X) = I 0 exp[-o(X)N£] 
such that the spatially averaged cross section is: 

tf(X) = ~^ln Jiexpf-Ojj (X)Nl] + ^exp[-a a (X)Nf]} (5-2) 

Integrating eq. 5-2 over the entire manifold gives the integrated cross section: 

f = ~ J lnjlexpt-a^Nf] + |exp[-a a (X)W]} d X (5-3) 

manifold 
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Figure 5-2a Visible Absorption Spectrum in YLF:Ho ic-polarization 




62 



Cross section (xlO ' 20 cm 2 ) Cross section (xlO 


3.0 



n 


2.5 - m 
2.0 - 

1.5 
1.0 
0.5 
0.0 

300 350 400 450 500 550 600 650 700 

Wavelength (nm) 



Figure 5-3a Visible Absorption Spectrum in YLF:Ho o-polarization 



Figure 5-3b Infrared Absorption Spectrum in YLF:Ho G-polarization 
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Equation 5-3 was used for the calculation of integrated absorption profiles. These 
were then substituted into eq. 2-13, along with the other constants and parameters 
appearing in that equation, to find the experimental electric dipole line strengths. 

A computer program was written to perform a least squares fit between the 
experimental linestrengths found using eq. 2-13, and the theoretical linestrengths of eq. 2- 
5. The object is to vary the Judd-Ofelt intensity parameters until the difference in the 
squares of the experimental and theoretical line strengths is minimized. The procedure is 
the same as that outlined by Krupke [Kr71], If the three Judd-Ofelt parameters are made 
to form a column vector ft , and the doubly reduced matrix elements form an n x m matrix 
M, where n is the number of transitions to fit, and m is the number of parameters to vary 
in the fit (in this case m=3 corresponding to the 3 Judd-Ofelt parameters), then the 
theoretical line strength then be written as a 1 x n column vector. This can be expressed 
as: 

S = MQ (5-4) 

The set of Judd-Ofelt parameters that minimizes the sum of the squares of the 
differences between the theoretical and experimental line strengths is found from matrix 
algebra to be (Appendix B): 

£i° = (5-5) 

where Mt is the adjoint of matrix M and (MtM)' 1 is the inverse matrix of MtM. The 
measure of the quality of the fit is given by the root mean square deviation, expressed as: 

RMS AS = [ E(AS) 2 /(# Transitions to fit - # of parameters to fit)] 1/2 
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This procedure was also carried out by hand to ensure that the computer program 


was giving correct results. The Judd-Qfelt parameters for YLF:Ho were found to be: 


Q 2 = 0.7483x10*20 cm 2 
Q 4 = 1.6168 xlO -20 cm 2 
Q 6 = 1-2269 xlO’ 20 cm 2 


This set of parameters were utilized to calculate the manifold to manifold 
spontaneous emission probabilities A(yJ,Y'J') using eq. 2-5 and eq. 2-12. From the 
discussion in chapter 2, it is recalled that J refers to the initial state and J’ refers to the 
final state. The radiative lifetime is found by summing over the inverse of the transition 
probabilities for all states below the one under consideration. This is expressed as: 


T rad 


1 

XA0|fJ ,yT) 


(5-6) 


The branching ratio for a given transition probability can be written as: 


o _ ACyJ.yT) 
Pji " IA(\jrJ,yT) 


(5-7) 


The results of the Judd-Ofelt analysis are presented in table 5-2. The first part of 
the table shows the manifolds used in the fit, their corresponding square matrix elements 
and the measured and theoretical values for the electric dipole linestrength. The 5 l 7 
manifold was excluded from the fit due to its strong magnetic dipole contribution. The 
second part of the table shows the values found for the transition probability (Aed) from 
the initial J state to the final J' state, the branching ratio for the given transition, and the 
radiative lifetime from the J' state to the ground state found from eq. 5-6. 
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Table 5-2 Judd-Ofelt Fit for YLF: Ho 


Transition 
from 518 

fU (2)]2 



X(nm) 

Smeas. 

xlO‘ 20 cm 2 

S theor. 
xl0' 20 cm 2 

■ 1 1 ! T 1 1 — 

■ ll" — 

■hi 



0.1421 j 

■l^«i 

3L9+5G3 

0.0185 

0.0052 


345 

0.2118 

0.2275 

3D2+3H6+5G5 

0.2155 

0.1969 


360 

0.9189 

0.6858 

3K7+5G4 

0.0058 

0.0361 


385 

0.17% 

0.1483 

3G5 

0 

0.5338 


418 

0.9251 

0.8620 

5F1+5G6 

1.5201 

0.8410 

0.1411 

448 

2.6378 

2.6704 

3K8 

0.0208 

0.0334 

0.1535 

467 

0.2814 

0.2583 

3F2 

0 

0 

0.2041 

473 

0.2789 

03510 

3F3 

0 

0 

03464 

486 

0.4729 

0.4260 

5F4+5S2 

0 

0.2392 

0.9339 

540 


15348 

5F5 

0 

0.4250 

05687 

645 


13856 

515 

0 

0.0100 

0.0936 

886 


0.1312 

516 

0.0084 

0.0386 

0.6921 

1175 

i 0.8090 

0.9198 


Judd-Ofelt Parameters: G 2 = 0-7483xl(r 20 G 4 = 1.6168xl(T 20 fl 6 = 1.2269xl0' 20 

Root Mean Square deviation: 8 = 0.10909xl0* 2 ® 
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5.3 Emission Cross Sections in Y I F ;Hq 


The methods of obtaining emission cross sections have already been discussed in 
chapter 2. The emission cross sections for YLF:Ho can be obtained directly from eq. 2- 
17. Since YLF is an anisotropic crystal, eq. 2-17 must be rewritten to account for the 
different polarizations of the measured luminescence spectrum. This is easily done by 
weighting the lineshape function as 2/3 o-pol. and 1/3 Jt-pol. Equation 2-16 for a uniaxial 
crystal becomes : 


g“(v) = 


3I a a) 

7 ? 


(5-8) 


where a denotes a given polarization (jc or a). 
Equation 2-17 now becomes: 


g“(\) = 


8jtcn 2 x rad 


_3IoO)__ 

2Iq(X)+I,(X) a 


(5-9) 


Emission spectra were taken for singly doped YLF:Ho. Since Holmium shows no 
absorption in the range accessible by our diode laser (785 nm), the crystal was pumped 
with a 35mW Helium-Neon laser operating at 635 nm. This wavelength corresponds to 
the 5 F 5 manifold in YLF:Ho (see figure 5-1). Emission spectra were taken over both 
polarizations in the wavelength range 800 - 2200 nm, and then corrected for grating 
efficiency and detector response with a tungsten lamp source operating at 2854°K. 
Calibrating the spectra is a very important procedure because over a given wavelength 
range the grating has a given efficiency curve for the scattering of light and the detector 
used has a response curve associated with it. As a result, the luminescence spectra 
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collected do not always reflect the relative intensities of the peaks in a wavelength range 
encompassing a manifold. This procedure seems to be one which is passed along to 
researchers through experience. I was not able to find the procedure in print anywhere. 
Therefore, for the benefit of researchers who are unaware of this technique, it is included 
in Appendix C. 

Emission cross sections were calculated from eq. 2-24 using absorption cross 
section data (see figures S-2 and 5-3), and calculated partition functions from energy level 
data [KaWo76], Table 5-3 contains Boltzmann fractions for the Stark levels of the four 
lowest manifolds in YLF:Ho. The partition functions for each manifold appear at the 
bottom of the tables. The emission cross section was determined independently from the 
modified Fuchtbauer-Ladenberg equation (eq. 5-9) by adjusting the radiative lifetime such 
that values for the emission cross section were in best agreement with those obtained 
from eq. 2-24. In this way a set of radiative lifetimes for the manifolds studied in YLF:Ho 
were determined. The results of this procedure are presented in figures 5-4, 5-5, and 5-6. 
The radiative lifetimes for the manifolds studied are given in table 5-4. The branching 
ratios for the appropriate manifolds were found from the Judd-Ofelt results in table 5-2. 
The radiative lifetime found for the manifold in YLF:Ho is in good agreement with the 

value 14.2 ms found by [Pa92], 

5.4 l umi n escenc e Pecay in YLEiHo 

A crystal of singly doped YLF.Ho with 0.5% Ho concentration was excited using 
a Ti:Al203 (Titanium Sapphire) laser operating at 25 mJ with 100 ns pulses. The 
wavelength of the laser was tuned to 892 nm, which corresponds to the 5 Is manifold of 
YLF:Ho (see figure 5-1). Emission was monitored with a Hewlett-Packard 54504 A 
digitizing oscilloscope at the wavelengths 537nm ( 5 S2), 658 nm ( 5 Fs), and 912 nm ( 5 Is). 
The data were collected with the Labview software package. 
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Table 5-3 Boltzmann Fractions in YLF :Ho @300K 


Level# 

Ho 518 
Energy 
(cm -1 ) 

Ho 518 
Boltz. 
Irac. 

Ho 517 
Energy 
(cm' 1 ) 

Ho 517 
Boltz. 
frac. 

Ho 516 
Energy 
(cm' 1 ) 

Ho 516 
Boltz. 
frac. 

Ho 515 
Energy 
(cm' 1 ) 

Ho 515 
Boltz. 
frac. 

1 

0 

0.10871 

5153 

0.08760 

8670 

0.09298 

11240 

0.10460 

2 

0 

0.10871 

5157 

0.08635 

8671 

0.09253 

11240 

0.i046Q 

3 

7 

0.10512 

5157 

0.08635 

8680 

0.08862 

11240 

0.10460 

4 

23 

0.09735 

5164 

0.08310 

8680 

088862 

11245 

0.10212 

5 

48 

0.08635 

5175 

0.07771 

8686 

0.08561 

11248 

0.10066 

6 

56 

0.08310 

5185 

0.07514 

8686 

0.08561 

11248 

0.10066 

7 

72 

0.07696 

5185 

0.07514 

8687 

0.08570 

11252 

0.09875 

8 

71 

0.07696 

5207 

0.06794 

8696 

0.08208 

11298 

0.07920 

9 

217 

0.03839 

5229 

0.06114 

8701 

0.08013 

11327 

0.06891 

10 

270 

0.02978 

5229 

0.06114 

8768 

0.05811 

11327 

0.06891 

11 

270 

0.02978 

5233 

0.05912 

8783 

0.05407 

11333 

0.06696 

12 

276 

0.02893 

5291 

0.04498 

8783 

0.05407 



13 

283 

0.02798 

5293 

0.04476 

8796 

0.05081 



14 

290 

0.02705 

5293 

0.04476 





15 

303 

0.02542 

5293 

0.04476 





16 

303 

0.02542 

j 






17 

315 

0.02400 







Part. 

func. 

Z= 9.1990 

Z= 11.4150 

Z= 10.7550 

Z= 9.5601 


Table 5-4 Radiative Lifetimes From Emission Cross Sections 


Transition 

Branching ratio § 

Trad/P 

Trad 

5 l5 -> 5 l8 

§58 = 0.408 

16.1 ms 

6.56 ms 

5 l6-> 5 l8 

§68 = 0.914 

6.20 ms 

5.66 ms 

00 
K— t 

* 

r- 

7? 

§78 = 1.000 

14.0 ms 

14.0 ms 


69 

















































































































































>ss section (xl 


7c-pol. absorption 
spectrum 


Derived rc-pol. 
emission spectrum 


o-poL absorption 
spectrum 



920 930 


Derived o-pol. 
emission spectrum 



900 910 920 930 


rc-pol. emission spectrum 

(c.) 

^*4- = 16.1 ms . 


Pss l 

yV_/ 

1 




o-pol. emission spectrum 

^ = 17.8 ms 

P58 


Wavelength (nm) 


Wavelength (nm) 


Figure 5-4 Cross Sections of the 5 Is Manifold in YLF:Ho 






























The decay curves are shown in figure 5-7. Figure 5-7a shows the decay from the 
5 I 5 manifold at 912 nm. This manifold showed only one exponential decay with a decay 
time of 20ps. Figure 5-7b shows the decay from the 5 Fs manifold at 658 nm. This 
manifold showed an exponential rise with a rise time of 27ps and an exponential decay 
with a decay time of 105ps. Figure 5-7c shows the decay from the 5 S2 manifold at 537 
nm. This manifold showed an exponential decay with a decay time of 105ps. From these 
data, luminescence lifetimes at room temperature (295K 0 ) can be assigned. The lifetime of 
the 5 l 5 manifold is 20ps. The lifetime of the s Fs manifold is 27ps and shows up as a rise 
time in Fig. 5-7b since this level if fed by a higher lying level, namely the 5 S2- The lifetime 
of the 5 S 2 manifold is 105|is which show up as a decay in both figures 5-7b and 5-7c. It 
should be pointed out that the mechanism by which the 5 S 2 manifold becomes populated 
when pumping the lower energy 5 Is manifold is not clear at this time. At any rate, these 
assignments of lifetimes are consistent with measurements conducted by exciting the 5 S 2 
manifold and the 5 Fs manifold directly using a tunable dye laser. Table 5-5 is a tabulation 
of these results. 


Table 5-5 Lifetimes Observed Under Different Excitation Wavelengths 



Emission from 

5 S2__ 

5F5 

5j 5 

5 S 2 @535 nm 

Rise time (ps): 


27.3 

32.0 


Decay time (ps): 


115 

100 

5 F 5 @635 nm 










5 l 5 @890 nm 
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Normalized Intensity 




Timers) 

Figure 5-7 Decay Curves in YLF:Ho Exciting the 5 Is Manifold @892 



Lifetime measurements have also been conducted on the same YLF:Ho sample 
using as an excitation source a CoMgF 2 (Cobalt Magnesium Fluoride) laser tuned to 
1.96pm to directly excite the 5 l 7 manifold of YLF:Ho. The emission from this manifold 
was monitored at 2.05p.m. The decay curve is shown in figure 5-8. The 5 l 7 manifold 
showed a single exponential decay with a decay time of 16.5 ms. 



Figure 5-8 Decay from YLF:Ho Exciting the 5 l 7 Manifold @ 1.96pm 
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CHAPTER 6 

YLF: Tm - Results and Interpretation 

We dance round in a ring and suppose, 
but the Secret sits in the middle and knows. 

Robert Frost 

The Secret Sits 


6.1 Energy Levels of YLF: Tm 


The energy levels for Tm 3+ in UYF4 have previously been measured and are 
reported in [Je71] and [Du91], The energy levels presented in these two references show 
some differences. Our low temperature absorption measurements made on YLF.Tm show 
better agreement with [Du91], so the experimental set of energy levels there will be used 
exclusively. An energy an energy level diagram of YLF:Tm is shown in fig. 6-1. 



25,000 

1— M 
• 

20,000 

a 

0 

15,000 

fa 


D 

G 

10,000 

m 


5000 


0 




Figure 6-1 Energy Level Diagram of YLF .Tm 
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A set of tables for the upper and lower transition wavelengths and energies for a 
given manifold were made from the energy levels reported in [Du91]. The transition 
wavelengths in nm are given in table 6-la , and for transition energy in cm* 1 in table 6-lb. 


Table 6- la Range of Transitions (nm) in YLF: Tm 



3H6 

3F4 

3H5 

3H4 

3F3 

3F2 

1G4 

1D2 

116 

3H6 

0 

1674 

1171 

775 

685 

654 

463 

356 

285 


0 

1930 

1271 

821 

709 

681 

486 

363 

291 

3F4 

- 

0 

3405 

1371 

1111 

1033 

626 

444 

340 



0 

4325 

1508 

1169 

1096 

666 

454 

347 

3H5 

- 

- 

0 

2170 

1584 

1430 

753 

505 

374 


- 

■- 

0 

2460 

1670 

1524 

803 

514 

381 

3H4 

- 

- 

- 

0 

5005 

3736 

1115 

646 

446 


- 

- 

- 

0 

6138 

4539 

1237 

663 

457 

3F3 

- 


- 

- 

0 

13245 

1420 

737 

488 


- 


- 

- 

0 

20120 

1568 

748 

496 

3F2 

- 


- 

- 

- 

0 

1546 

770 

502 


- 


- 

- 

- 

0 

1755 

788 

514 

1G4 

- 


- 

- 

* 

- 

0 

1408 

713 


- 

- 

- 

- 

- 

- 

0 

1562 

759 

1D2 

- 

■*- 

- 

- 

- 

- 

- 

0 

1421 


- 

- 

- 


- 

- 

- 

0 

1502 

116 

- 


- 

- 

- 

- 

- 

- 

0 



Table 6-lb. Range of Transitions (cm- 1 ) 

of YLF: Tm 

0 


.2M 2E4 2H5 3H4 3F3 3F2 1G4 1 B2 ll£ 


3H6 

0 

5972 

8535 

12891 

14597 

15275 

21562 

28075 

34998 


0 

5180 

7865 

12180 

14101 

14675 

20554 

27542 

34310 

3F4 

- 

0 

2936 

7292 

8998 

9676 

15963 

22476 

29399 


- 

0 

2312 

6627 

8548 

9122 

15001 

21989 

28757 

3H5 

- 

- 

0 

4607 

6313 

6991 

13278 

19791 

26714 


- 

- 

0 

4064 

5985 

6559 

12438 

19426 

26194 

3H4 

- 

“ 

- 

0 

1998 

2676 

8963 

15476 

22399 


- 

- 

- 

0 

1629 

2203 

8082 

15070 

21838 

3F3 

- 

- 

~ 

- 

0 

755 

7042 

13555 

20478 


- 

- 

- 

- 

0 

497 

6376 

13364 

20132 

3F2 

- 

- 


- 

- 

0 

6468 

12981 

19904 


- 

- 

- 

- 

- 

0 

5698 

12686 

19454 

1G4 

- 

- 


- 

- 


0 

7102 

14025 


- 

- 


- 

- 

— 

0 

6399 

13167 

1D2 

- 

- 

- 

- 

- 

- 

• 

0 

7037 


- 


*- 

- 


- 

- 

0 

6654 

116 

- 



- 

- 

- 

- 

- 

0 


- 

- 



- 

- 

- 

- 

0 
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6.2 ludd-Ofelt analysis of YLF;Tm 


A theoretical discussion of the Judd-Ofelt theory has been given in chapter 2, and 
its application was described in chapter 5 for YLF:Ho. The result of the Judd-Ofelt 
analysis for YLF:Tm is given here. 

Absorption measurements were made on a sample of Tm YLF doped with 0.5% 
thulium. Spectra were taken for both polarizations in the wavelength range 300 - 2300 
nm. The absorption spectra for jc-polarization appear in figures 6-2, and the absorption 
spectra for o- polarization appear in figures 6-3. 

As was discussed in the last chapter, equation 5-3 was used for the calculation of 
integrated absorption profiles, the values of which were then substituted into eq. 2-13 to 
find the experimental electric dipole linestrengths. 

The method of least squares was applied in the same way as was done for 
YLF:Ho in chapter 5. The Judd-Ofelt intensity parameters were found for YLF:Tm to be: 

Q 2 = 1 6673 xlO* 20 cm 2 
Q 4 = 1 .9184 xlO* 20 cm 2 
Q 6 = 1.0897 xlO' 20 cm 2 

These intensity parameters were used to calculate the manifold to manifold 
spontaneous emission probabilities. The results are presented in table 6-2. The meaning 
of this table is the same as discussed in chapter 5. The 3 Hs was excluded from the fit due 
to its strong magnetic dipole contribution. It should be expected that the error associated 
with the spontaneous emission probabilities for YLF:Tm will be larger than that for 
YLF: Ho. There are a couple of reasons why this is believed to be true. The first concerns 
the width of the manifolds in Tm (e.g., the energy difference between the highest stark 
level and lowest stark level). The use of eq. 2-12 involves the use of an average 
wavelength for a given manifold. Since the widths of the manifolds are about twice as large 
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as those in Holmium, then there could possibly be more error associated with the 
calculations. The second reason has to do with the number of manifolds available for the 
least squares fit. As can be seen from table 6-2 there were only five manifolds available 
for a fit in which three parameters were being varied. This is a disadvantage since it 
requires more accuracy in the measured values to obtain a good set of fitting parameters. 
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Cross section (xlO cm 2 ) Cross section (xlO 



Figure 6-2a Visible Absorption Spectrum in YLF:Tm Jt-polarization 



Wavelength (nm) 


Figure 6-2b Infrared Absorption Spectrum in YLF:Tm 7t-polarization 
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Cross section (xl0‘ cm 2 ) Cross section (xlO 



Figure 6-3a Visible Absorption Spectrum in YLF:Tm o-polarization 



Figure 6-3b Infrared Absorption Spectrum in YLF:Tm o-polarization 
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Table 6-2 Judd-Ofelt Fit for YLF:Tm 


Transition 
from 518 

[U<2>]2 

[U(4)]2 

[U<6>]2 

X(nm) 

5 meas. 
xl0" 20 an 2 

S theor. 
xl0‘ 20 cm 2 

*D2 

b 

05156 

0.0928 

357 

0.7209 

0.7065 

1G4 

0.0483 

0.0748 

0.0125 

469 


0.2371 

3F2+3F3 

0 

0.3164 

1.0992 

672 


1.8047 

3H4 

0.2373 

0.1090 

05947 

784 

1.2633 

1.2504 

3F4 

05375 

0.7261 

02382 

1725 

25347 

2.5433 


Judd-Ofelt Parameters: fi 2 = 1. 6673x10-20 fl 4 = 1.9184xl(r 20 ft 6 =1-0897x10-20 

Root Mean Square deviation: 8 = .026938x10*20 


Transition 

J->F 

wavelength 

Xfnm) 

S theor. 
xlO' 20 cm 2 

Aed 

(sec -1 ) 

Amd 

(see -1 ) 

branching ratio 

P 

Trad <Fs) 

1D2-1G4 

1486.3 

0.675 

79.18 


0.006 


1D2-3F2 

7885 

0.693 

587.65 


0.044 


1D2-3F3 

745.2 

0.401 

38753 


0.029 


1D2 - 3H4 

657.0 

0.473 

669.11 


0550 

75 

1D2-3H5 

504.1 

0.023 

71.49 


0.005 


1D2-3F4 

450.4 

1.143 

5092.22 


0584 


1D2-3H6 

357.1 

0.702 

6373.78 


0.481 


1G4-3F2 

1633.7 

0.211 

10.63 


0.007 


1G4-3F3 

14945 

0528 

34.75 


0.022 


1G4-3H4 

1177.4 

0.113 

1527 


0.010 


1G4-3H5 

7625 

0.755 

37759 


0244 

645 

1G4-3F4 

6462 

0.729 

602.71 

/ 

0589 


1G4-3H6 

470.1 

0235 

50951 


0529 


3F2-3F3 

17513.1 

0.151 

0.02 


0.000 


3F2-3H4 

4215.9 

0.869 

450 


0.006 


3F2-3H5 

1431.0 

1.192 

160.91 


0207 

1289 

3F2-3F4 ' 

1069.1 

0.671 

21825 


0281 


3F2-3H6 

660.0 

0.281 

39229 


0506 


3F3-3H4 

55525 

1.082 

1.49 


0.001 


3F3-3H5 

1558.4 

1.708 

12759 


0.084 


3F3 - 3F4 

1138.6 

0.185 

35.61 


0.023 

660 

3F3-3H6 

685.9 

1532 

1351.08 


0591 


3H4-3H5 

2166.4 

0.943 

2021 


0.030 


3H4-3F4 

14325 

0.701 

52.46 


0.078 

1480 

3H4-3H6 

8825 

1502 

603.06 


0592 


3H5 - 3F4 

42265 

1573 

3.06 


0523 


3H5-3H6 

1225.0 

1516 

12955 

46.0 

0.977 

5614 

3F4-3H6 

1725.0 

2.499 

106.66 


1500 

9375 
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6.3 Emission Cross Sections in YLF;Tm 


Emission spectra were taken for singly doped YLFiTm. Thulium shows strong 
absorption in the range accessible by our diode lasers. The crystal was pumped with a 
300mW diode laser operating at 785 nm. This wavelength corresponds to the 3 H 4 
manifold in YLFiTm (see figures 5-2 and 5-3). Emission spectra were taken over both 
polarizations in the wavelength range 400 - 2200 nm, and then corrected for grating 
efficiency and detector response with a tungsten lamp source operating at 2854 °K. 

The emission cross sections for YLFiTm were calculated in the same way as was 
done for YLFiHo (see chapter 5, section 3 for details). The energy level data needed for 
the calculation of the partition functions was found from [Du91], Table 6-3 contains 
Boltzmann fractions for the Stark levels of the four lowest manifolds in YLF:Tm. The 
partition functions for each manifold appear at the bottom of the table. No fluorescence 
was detected from the 3 Hs manifold, so only the results for the 3 H 4 and 3 F 4 manifolds are 
presented in figures 6-4 and 6-5. 


Table 6-3 Boltzmann Fractions in YLF:Tm @300K 


Level# 

Tm3H6 

Energy 

(cm -3 ) 

Tm3H6 

Boltz. 

frac. 

Tm3F4 

Energy 

(cm -3 ) 

Tm3F4 

Boltz. 

frac. 

Tm3H5 

Energy 

(cm -3 ) 

Tm3H5 

Boltz. 

frac. 

Tm3H4 

Energy 

(cm -3 ) 

Tm3H4 

Boltz. 

frac. 

1 

0 

0.19203 

5599 

0.29133 

8284 

0.16844 

12599 

0.18981 

2 

30 

0.16630 

5757 

0.13655 

8300 

0.15600 

12624 

0.16836 

3 

30 

0.16630 

5757 

0.13655 

8300 

0.15600 

12643 

0.15370 

4 

56 

0.14680 

5765 

0.13141 

8319 

0.14241 

12643 

0.15370 

5 

270 

0.05260 

5820 

0.10094 

8501 

0.05949 

12745 

0.09424 

6 

305 

0.04447 

5942 

0.05623 

8501 

0.05949 

12804 

0.07101 

7 

319 

0.04158 

5968 

0.04963 

8519 

0.05457 

12835 

0.06120 

8 

334 

0.03870 

5972 

0.04869 

8531 

0.05152 

12835 

0.06120 

9 

334 

0.03870 

5972 

0.04869 

8534 

0.05078 

12891 

0.04678 

10 

372 

0.03225 



8534 

0.05078 



11 

407 

0.02727 



8535 

0.05054 



12 

407 

0.02727 







13 

419 

0.02574 







fgQM 

ES9H 

Z= 5.2074 

Z= 3.4326 

Z= 5.9369 

Z= 5.2684 
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Figure 6-4 Cross Sections of the 3 Ht Manifold in YLF:Tm 



















6.4 Luminescence Decay in YLF:Tm 


Crystals of singly doped YLF.Tm were excited using a Ti:Al203 (Titanium 
Sapphire) laser tuned to 780 nm, corresponding to the 3 H4 manifold (see figure 6-1). 
Emission was monitored with a Hewlett-Packard 54504A digitizing oscilloscope at 812 
nm. The data was collected with the Labview software package. Emission was also 
monitored with a LeCroy CAMAC at 1.7pm for excitation at 780nm and 1.7p.m. This 
data was collected with the LeCroy Catalyst software package. The details of the 
experimental apparatus used for these experiments was given in chapter 4. 

The decay of Tm ions from the 3 H4 manifold exhibit a non-exponential decay 
behavior that strongly shortens with increasing Tm concentration. The shortening of the 
decay with increasing Tm concentration is a result of quenching via cross-relaxation. In 
this process an ion excited into the 3 H4 manifold interacts with a nearby ground state ion 
in 3 H6 and undergoes an exchange of energy, resulting in two ions in the 3 F4 manifold. 
The cross-relaxation mechanism, interpreted as the transfer Tm( 3 H4 -> 3 F4);Tm( 3 H6 -> 
3 F4 ), is a phonon assisted transfer in the case of YLF since the gaps 3 H4 - 3 F4 = 5972 
cm' 1 and 3 F4 - 3 H6 = 6627 cm* 1 exhibit no overlap. 

Measurements of the decay from the 3 H4 manifold of Tm after excitation of the 
3 H4 with a Ti:Al2C>3 laser tuned to 780 nm were made on samples of Tm(.5) YLF, Tm(3) 
YLF, and Tm(12) YLF. The Tm(.5) YLF sample exhibited exponential decay with a 
decay time of 2000ps. The Tm(3) YLF and Tm(12) YLF samples were non-exponential 
in behavior and exhibited increasingly faster decay with increasing Tm concentration. The 
shortening of the decay is a result of the cross-relaxation process just discussed. In 
considering this process, we must also consider that the cross-relaxation can be enhanced 
by migration among the Tm 3 H4 ions described by the process Tm( 3 H4 -> 3 H6);Tm( 3 H6 
-> 3 H4 ). The migration is usually treated as a diffusion process or a hopping process. In 
the diffusion process, the excitation makes many jumps between donors before being 
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deactivated by an acceptor. On the other hand, in the hopping process, the excitation 
makes a single jump in the vicinity of an acceptor and is subsequently deactivated by the 
acceptor. The diffusion approach is warranted if the average transfer time between donors 
is much greater than the transfer time between the donor and acceptor, whereas the 
hopping approach is warranted if the donor-donor transfer time is less than or equal to 
the donor-acceptor transfer time. These statements can be quantified in terms of the 
microscopic interaction parameters as: Cdd « Cda for diffusion migration and Cdd - 
Cda for hopping migration [608 1], In general, at low concentrations of active ions in the 
crystal the concept of defining a diffusion constant becomes problematic because of the 
variation of ion-ion transfer rates in different environments. In other words, the problem 
becomes more discrete in nature and the diffusion gradient cannot be defined. In such 
cases, the hopping model can be used. 

The dopant ion concentrations used in this study are not extremely large, so the 
hopping mechanism is a likely choice for the interpretation of migration enhanced cross- 
relaxation. Before proceeding with this line of thought, however, there is a simple way to 
check whether the hopping interpretation is applicable in this case. The microscopic 
interaction parameters for the D->A and D->D transfer can be calculated easily using 
Dexter's theory [DeS3], It should pointed out that the donors and acceptors in this case 
are the same ion, namely Tm, but we will still write the interaction parameters as Cda to 
refer to cross-relaxation and Cdd to refer to migration. 

The interaction parameter for energy transfer in Dexter's theory is: 



(6-1) 
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where X=A for cross-relaxation and X=D for migration. The other terms in this equation 
have the usual meaning. 

Calculating the overlap integral and integrated absorption cross section in eq. 6-1, 
we found the interaction parameters for the Tm-Tm cross-relaxation and migration to be: 

C D A = 3.7083x10-47 C m 6 /ps 
C D d = 1.2624x10-45 cm 6 /ps 

Clearly, these values satisfy the relation Cdd ^ Cda- So, our reasoning for the use 
of the hopping mechanism over diffusion is validated. We can now proceed to interpret 
the luminescence decay from the Tm 3 H4 in terms of a hopping migration assisted cross- 
relaxation model. 

Burshtein [Bu72] is usually attributed with the initial development of the hopping 
mechanism of energy transfer, however, Artamonova's paper [Ar72] appeared in the same 
journal as Burstein's paper, but two months earlier. The paper by Artamanova is more 
lucid, so we will refer to his paper in the discussion that follows. 

In the hopping model, the time dependence of the donor decay can be described 
by the equation: 

|E« =exp [_X_V / 2. Wt ] (6-2) 

where to is the donor lifetime (radiative + nonradiative), y is a term describing the Tm-Tm 
cross-relaxation, and W is a term describing the migration enhanced cross relaxation. The 
y term for dipole-dipole interaction between ions is the same term that appears in the 
Forster-Dexter, Inokuti -Hi ray ama theories [De53, InHi65], It is given by: 
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(6-3) 


Y = 4k 3/2 N a <C da ) 1/2 

The term describing the migration enhancement of the cross-relaxation was first 
derived in [Bu72, Ar72j. Under the assumption of a dipole-dipole interaction in the 
hopping mechanism this term is: 

w = k(^) 5/2 N a N d (C da ) 1/2 (C dd ) 1/2 (6-4) 

The addition of a term which takes into account the migration among donors 
provides a way for the excitations in the Tm 3 H 4 to redistribute themselves randomly at 
some hopping rate Th' 1 - This causes an enhancement of the cross-relaxation because 
excitations which are relatively isolated can migrate until they find a quenching center. 

Equation 6-1 was fit to the experimental decay curves of the 3 H 4 in samples of 
Tm(3)YLF and Tm(12)YLF. The results appear in table 6-4 along with the values found 
from Dexter’s theory for comparison. 


Table 6-4 Values of Cda and Cdd in Tm YLF 



Cda (cm 6 /ns) 

CDD(cm 6 /ps) 

Cdd/Cda 

Tm(3)YLF 

(3.7591±.0924)xl0* 47 

(2.9091±.0100)xl0* 45 

0.0129 

Tm(12)YLF 

(6.7532±.0268)xl0- 47 

(4.0010±.0193)xl0 -45 

0.0168 

Dexter Theory 

3.7083X10* 47 

1.2624x10-45 

0.0293 


The fluorescence quantum efficiency can be found from the experimental decay 
curves by: 

Tlf = ^ J~[n D (t) / n D (0)]dt (6-5) 
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where x r is the radiative lifetime, taken to be 2400}is as found in section 6.3. Caird, et al. 
[Ca91] has integrated eq. 6-5 using eq. 6-2 as the integrand. The fluorescence quantum 
efficiency was found to be: 

% = Jt 1/2 xexp(x 2 )[l- erf(x)]} (6-6) 

where x r is the radiative lifetime, Xo is the lifetime in the absence of cross-relaxation 
(radiative + nonradiative), and x = y / 2-^(1 / T 0 + W) . The values for y and W can be 

calculated from equations 6-3 and 6-4 using the values of the interaction parameters found 
in table 6-4. The value for Xo is taken to be 2000ps, the decay time in the low 
concentration Tm(.5)YLF sample. Based on the values calculated from the quantum 
efficiency, we can estimate the cross-relaxation time, since the quantum efficiency reflects 
the ratio of the mean decay time and the radiative time of decay . That is, we can write: 


*cr 


= 'Hf 
“V ^0 


(6-7) 


A tabulation of the results appears in table 6-5. Unfortunately, we do not have 
enough singly doped YLF:Tm samples to make an assessment of the concentration 
dependence of the cross-relaxation. 


Table 6-5 Fluorescence Quantum Efficiencies in Tm YLF 



y(ps- 1/2 ) 

mm 

m ' 

Xcr(lts) 


0.0188 

0.0011 

0.2388 

n 

mm 

0.0101 

0.0284 

0.0114 

ISEfll 
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In contrast to YLF:Tm, crystals of YAG:Tm do not seem to exhibit migration 
enhancement of the cross-relaxation. The decay of the Tm 3 H 4 manifold in YAG is well 
described by eq. 6-2 with W = 0. This is the case of static cross-relaxation described in 
the Forster-Dexter, Inokuti-Hirayama theories. The interaction parameter for YAG:Tm 
was found to be on the order of 10' 45 cm 6 /fis. This is two orders of magnitude greater 
than in YLF:Tm, indicating that the interaction between ions for cross-relaxation is much 
stronger in YAG than in YLF. This is supported by observation of the decay curves of 
the Tm 3 H 4 manifold in YAG and YLF. The lifetime in YAG quenches much faster with 
increasing Tm concentration than YLF. It seems to be the case, then, that migration among 
Tm ions can be neglected in YAG, but must be considered in YLF in order to correctly 
model the decay behavior. 
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CHAPTER 7 

YLF: Tm,Ho - Results and Interpretation 


Nature uses only the longest threads to weave her patterns, 
so that each small piece of her fabric reveals the organization 
of the entire tapestry. 

Richard Feynman 

The Character of Physical Law 

7.1 Sensitised luminescence in YLF-i.TmJfo 

In the previous two chapters the focus centered on the spectroscopy and 
excitation dynamics of YLF crystals singly doped with either Thulium ions or Holmium 
ions. This chapter will concentrate on YLF systems co-doped with both Thulium and 
Holmium ions. Sensitized luminescence involves the initial excitation of a set of ions, 
termed sensitizers or donors (D), followed by the transfer of energy to other ions termed 
activators or acceptors (A), and the subsequent emission of light by acceptor ions. This 
description is very general, but represents many possible processes responsible for the 
transfer of energy from donor to acceptor. These include resonant transfer, non-resonant 
transfer, cross relaxation and upconversion. The 2.0pm laser emission from diode pumped 
Thulium sensitized Holmium has been the subject of numerous studies for various crystal 
hosts. Two of the most popular host materials studied have been Y 3 AI 5 O 12 (YAG) 
(Fa87, Fa88, Km94], and LiYF 4 (YLF) [He89, Es90, Bu92, St93], It is, therefore, useful 
to examine some of the spectroscopic features of diode pumped YLF:Tm,Ho. 

YLF:Tm,Ho crystals were excited with a diode laser operating at 785 nm, 
corresponding to the 3 H 4 manifold in Thulium. Emission from the crystals was monitored 
from 300-2300 nm. The transitions that were observed under these excitation conditions 
is pictured in Figure 7-1. The striking feature of this figure is the plethora of transitions 
that result from pumping one manifold in Thulium. Based on the values for the range of 
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Energy (x 10 c 


transitions of YLF:Ho in Table 5-la and YLF:Tm in Table 6-la, the following 
assignments of the numbered transitions in figure 7-1 were made. Figure 7-2 shows the 
scaled emission spectra depicting the relative strengths of these transitions. 



Tm Ho 


Figure 7-1 Observed Transitions in YLF:Tm,Ho Exciting the Tm 3 H 4 @785 nm 
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Scaled Intensity (A.U.) Scaled density (A.U.) 




la. ) Tm: *G4 -> 3 F4 (626-666 nm) 

lb. ) Tm: ^4 -> 3 H<> (463-486 nm) 

2.) Tm: 3 F 3 -> 3 H6 (684-708 nm) 

3a.) Tm: 3 H4 -> 3 H 5 (2170-2460 nm) 

3b.) Tm: 3 H4 -> 3 F 4 (1371-1508 nm) 

3c.) Tm: % -> 3 H6 (779-818 nm) 

4.) Tm: 3 F 4 -> 3 He (1889-2066 nm) 

5a.) Ho: 5 S 2 -> 5 l7 (730-750 nm) 

5b.) Ho: 5 S 2 -> 5 h (539-550 nm) 

6. ) Ho: 5 F 5 -> 5 l8 (638-659 nm) 

7. ) Ho: 5 I 5 -> 5 l8 (882-915 nm) 

8. ) Ho: 5 I 6 -> 5 l8 (1 136-1 196 nm) 

9. ) Ho: 5 I 7 -> 5 Is (1693-1924 nm) 

It is obvious from these observations that the YLF:Tm,Ho system presents a 
formidable problem toward understanding the modes of excitation and interaction of the 
various manifolds. The complexity of some of these processes precludes a rigorous 
analysis of the available experimental data. When this is the case, a more prudent 
approach to the interpretation of the data is to resort to a descriptive study of the 
phenomena. No attempt will be made here to entirely explain all of the observed features, 
but given a set of suitably designed experiments some insight can be gained toward this 
goal. It will be helpful to make some descriptive observations concerning Figure 7-1. 

The transitions pictured in figure 7-1 for YLF:Tm,Ho were the only ones 
observed in the range 300-2300 nm. The only two manifolds that were found to show no 
luminescence were the Ho 5 Lj and the Tm 3 Hs. In the case of the Ho 5 l4 this can be 
explained by its very small electric dipole transition probability (see table 5-2) and the 
small energy gap to the next lower lying manifold. Its primary channel of de-excitation is 
nonradiative and no optical radiation would be expected. In the case of the 3 Hs the most 
likely explanation is the following. The energy gap from the 3 H4 to the 3 Hs is large, 
precluding any substantial excitation of the 3 Hs from nonradiative decay of the 3 H4. The 
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feeding of the 3 Hs from the 3 H4 through optical means accounts for only 3% of the total 
3 H4 luminescence, as can be seen from the branching ratio in table 6-2. The low 
population in the 3 Hs in conjunction with a small energy gap to the next lower lying 
manifold, and a moderately large radiative lifetime would lead to very little luminescence. 
This situation is exacerbated further at higher Tm concentrations where de-excitation of 
the 3 H4 manifold is dominated by cross relaxation to 3 F4. This explanation assumes that 
the 3 H5 is not substantially populated by any energy transfer processes. 

One striking feature about Figure 7-1 is that more than half of the observed 
transitions originate above the pump manifold. This is indicative of upconversion 
processes. This conclusion is supported by observing the change in intensity of the 
transitions in these manifolds as a function of pump power. Exciting the 3 HL| manifold 
with a diode laser operating at 785 nm, the intensity of the Ho 5 S2 @550 nm. Ho 5 Is 
@890 nm. Ho 5 l7 @2.06 pm, and Tm *G4 @ 485 nm, Tm 3 H4 @812 nm, Tm 3 F4 @1.7 
pm were measured as a function of absorbed power by varying the laser diode pump 
power from 100% to 5%. The results of these measurements for a Tm(4)Ho(.5) YLF 
crystal are shown in figure 7-3. 

Three types of behavior were seen. The intensity followed either linear, super 
linear, or sublinear change with absorbed power. Referring to figure 7-3, it can be seen that 
the Ho 5 S2, Ho 5 Is, and Tm *G4 exhibit super linear behavior, the Tm 3 H4 exhibits linear 
behavior, and the Ho 5 l7, Tm 3 F4 exhibit sublinear behavior. The Tm 3 H4 behavior is as 
expected since the pump band reflects only the linear change in intensity with pump 
power. The change in intensity of the Ho 5 S2, Ho 5 l5, and Tm *G4 with absorbed power 
reflect a gain in the population as a result of upconverted energy into the manifold. The 
Ho 5 l7 and Tm 3 F4 reflect a loss in the population with absorbed power as a result of a 
depletion of energy, suggesting that these manifolds give up energy to support direct 
energy transfer processes and upconversion processes. 
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Figure 7-3 Luminescence Intensity vs. Absorbed Power in YLF:Tm(4)Ho(.5) 






The task of identifying the channels of these energy transfer processes is a 
difficult one which is far from being solved. Ideally, the best method of approach would 
be to excite many different manifolds and observe the luminescence and decay 
characteristics of the various other manifolds under different pumping conditions and over 
a range of temperatures. Such experiments should also be done with a suitable range of 
dopant ion concentrations. Only in this way can a comprehensive and thorough analysis 
be done that leads to solid conclusions. In practice, this is not always possible due to 
constraints on the availability of the necessary apparatus. Therefore, a reliance on studies 
available in the literature in conjunction with ones own studies must serve as a basis for 
any conclusions that are made. 

A number of emission spectra under diode pumping of the Tm 3 H 4 of 
YLF:Tm,Ho crystals with various dopant ion concentrations and identical path lengths 
along the pump direction have been examined. Based on these emission spectra, the 
following observations can be made concerning the change of the intensity of the emission 
with dopant ion concentration: 

1 . ) Ho 5 S 2 intensity decreases with increasing Ho concentration 

and decreases with increasing Tm concentration 

2. ) Ho s Fs intensity increases with increasing Ho concentration 

and increases with increasing Tm concentration 

3. ) Ho s l 5 intensity decreases with increasing Ho concentration 

and increases with increasing Tm concentration 

4. ) Ho 5 I^ intensity decreases with increasing Ho concentration. 

5. ) Ho 5 l 7 intensity increases with increasing Ho concentration. 

6) Tm *G 4 intensity increases with increasing Ho concentration. 

7. ) Tm 3 F 2 , 3 F 3 intensity decreases with increasing Ho concentration 

and decreases with increasing Tm concentration 

8. ) Tm 3 H 4 intensity decreases with increasing Tm concentration 

and is independent of Ho concentration. 

9. ) Tm 3 F 4 intensity decreases with increasing Ho concentration 

and increases with increasing Tm concentration. 


100 



These observations are general statements on the observed behavior of the 
intensity of the luminescence of various manifolds in the YLF:Tm,Ho system under 
excitation of the 3 H4 manifold of Tm. The actual behavior may be more subtle and depend 
on the ratio of the concentrations. At any rate we can attempt to draw some conclusions, 
based on these observations, regarding how these manifolds become populated. 

Exciting the 3 H4 manifold of Tm produces luminescence from the this manifold 
that decreases with increasing Tm concentration. This is explained by the interaction of a 
Tm ion in an excited 3 H4 state with a ground state ion in the 3 H 6 state, resulting in an 
exchange of energy that produces two ions in the 3 F4 state. This process is known as 
cross-relaxation and can be referred to by the notation Tm( 3 H4-> 3 F4);Tm( 3 H6- > 3 F4). 
Observations eight and nine corroborate this interpretation because as the Tm 
concentration is increased the interaction among the Tm ions increases. The cross- 
relaxation process works towards depleting the 3 H4 in favor of populating the 3 F4. Thus, 
as Tm concentration is increased, the 3 H4 intensity decreases while the 3 F4 intensity 
increases. These are exactly observations eight and nine. This is valid for both singly 
doped YLF:Tm and YLF:Tm,Ho. With the addition of Holmium, however, energy 
transfer from the Thulium to the Holmium ions becomes possible. In observation nine we 
noted that the intensity of the Thulium luminescence from 3 F4 decreased with Holmium 
concentration, indicating that ions in the 3 F4 interact with Ho ions and exchange energy, 
depleting 3 F4 states and populating Ho states. 

One such channel of energy transfer involves direct transfer from Tm 3 F4 to the 
Ho Tm( 3 F4-> 3 H6);Ho( 5 l8-> 5 l7), in which an excited 3 F4 ion drops to the 3 H 6 ground 
state and a Ho ion is excited from the 5 Is ground state to the excited 5 l7 state. This is 
supported by the decay kinetics of these two manifolds. When exciting the 3 F4 manifold, 
this manifold shows an initial fast decay (tens of microseconds), followed by a long decay 
(thousands of microseconds) when looking at the 3 F4 emission. When looking at the 5 H 
emission, an initial fast rise is seen, followed by a long decay equal to the long decay of 
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the 3 F4 emission. This behavior will be explained in more detail later. At this point we 
are only concerned with the qualitative behavior. The Tm 3 F*->Ho 5 l7 energy transfer is 
the first step from which all other Ho manifolds become populated. 

For instance, Fan et al. [Fa 88 ] has postulated that the luminescence in the Ho 5 Is 
results from the upconversion process Tm( 3 F4-> 3 H6);Ho( 5 l7 -> 5 Is). This seems to be in 
agreement with observations three and nine. Since the 3 F4 emission increases with Tm 
concentration then the 5 Is emission would be expected to increase also with increasing Tm 
since there are more excited ions in the 3 F4 state available for upconversion. Also, if this 
upconversion is present then a sublinear rise in intensity of the 5 l7 emission with 
absorbed power and a super linear rise in intensity of the 5 Is emission with absorbed 
power would be expected. This is, in fact, the case as evidenced by figure 7 - 2 . 

Tyminski, et al. [Ty 89 ] has attributed the 5 S 2 luminescence as resulting from the 
upconversion process Tm( 3 H4-> 3 H6);Ho( 5 l7 -> 5 S2). Unlike the 5 Is luminescence, the 5 S2 
emission decreases with increasing Tm concentration (observation one). Since the 3 H4 
emission decreases with Tm concentration, then the 5 S2 emission would be expected to 
also decrease with increasing Tm because there are less excited ions in the 3 H4 state (due 
to cross relaxation to the 3 F4) available for upconversion. Like the 5 Is , the 5 S2 intensity 
shows a super linear rise with absorbed pump power. 

There has been some controversy regarding the importance of these upconversion 
processes with regards to depletion of the metastable 5 l7 laser level. A depletion of the 
population in the laser level decreases the 2 \i laser output by increasing the threshold (the 
condition regarding the required population inversion for lasing to occur) [Sv 82 ]. It has 
generally been assumed that the 5 S2 upconversion is unimportant [Fa 87 , Ki 87 ], but no 
justification for this assumption has been presented. Figures 7-2 show the relative 
intensities of the emission from a sample of Tm( 4 )Ho(. 5 )YLF by exciting the 3 H4 
manifold with a diode laser. As can be seen from these figures, the intensity of the 5 S2 

emission is 100 times smaller than the 5 l7 emission, while the 5 Is emission is only about 
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10 times smaller. A possible explanation for this could be that the 5 S 2 does not have as 
large a population as the 5 Is. This would seem to indicate that the 5 S 2 upconversion is 
less significant than the 5 Is upconversion in depleting the laser level. On the other 
hand, the nonradiative decay could also play a part in deactivating the 5 S 2 - 

7.2 Luminescence Pccay in 

Many measurements have been made of the temporal decay in YLF:Tm,Ho with 
different Tm and Ho dopant concentrations under different excitation wavelengths. These 
measurements are very useful for identifying the channels of energy transfer and 
determining the rates at which they occur. For instance, measuring the decay time of a 
particular manifold in singly doped YLF:Tm and YLFiHo and then measuring the decay 
of the same manifold in the co-doped sample can indicate the presence of energy transfer 
or trapping of energy by observing the shortening or lengthening of the decay time. 

Two laser excitation sources were used to measure the temporal behavior in co- 
doped YLF:Tm,Ho. The Ti:Al 2 C >3 was tuned to 785nm exciting the 3 H 4 manifold of Tm 
and tuned to 890nm exciting the sis manifold of Ho. The CoMgF 2 laser was tuned to 
1.75jim to excite the 3 F 4 manifold of Tm and tuned to 1 .96pm to excite the 5 l 7 manifold 
of Ho. Using these two excitation sources we were able to selectively excite the above 
indicated manifolds and observe the decay of the fluorescence emitted by ions in various 
manifolds. The wavelength selection of the fluorescence was achieved with the use of 
narrow band filters when exciting with the CoMgF 2 laser, and with a spectrometer when 
exciting with the Ti:Al 2 C >3 laser. The reason for this is that the energy of the CoMgF 2 
beam was 6mJ while the beam energy of the Ti:Al 2 C >3 was 25mj. With the low energy of 
the CoMgF 2 laser it was necessary to use filters to get sufficient signal. The use of filters 
has the drawback that while much more signal is observed there is always some leakage of 
scattered light from the pump beam which can mix with the fluorescence signal. This is 
undesirable so all attempts were made to keep this to a minimum. No preamplifiers were 
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used to minimize signal distortion and a 50Q resistor was used to terminate the signal to 
the CAMAC to eradicate the RC time constant. Before proceeding, it should also be 
mentioned that the pulse length of the CoMgF2 laser was measured to be 480ns and the 
pulse length of the Ti:Al2C>3 laser was <100ns. 

The remainder of this chapter will be concerned with analyzing some of the energy 
transfer processes in detail by modeling the decay behavior. 

7.2.1 Cross-Relaxation Ho( 5 l5-> 5 I 7 );Tm( 3 H 6 -> 3 F 4 ) in YLF:Tm,Ho 

The decay of the 5 Is manifold in singly doped YLF:Ho has already been presented 
in chapter 5. Figure 5-7 shows the decay of the Ho 5 Is , Ho 5 Fs , and the Ho 5 S2 resulting 
from excitation of the Ho 5 Is at 890nm. From these measurements, the natural decay time 
of the Ho 5 l5 was found to be 20ps. In systems co-doped with Tm, however, the decay 
becomes non-exponential and the lifetime is found to be shortened. It seems evident that 
there is an interaction between Tm ions and ions in the excited Ho 5 Is manifold. The most 
likely process is the cross relaxation process Ho( 5 l5-> 5 l7);Tm( 3 H6 -> 3 F 4 ), which is a 
reverse process to the upconversion Tm( 3 F 4 -> 3 H6);Ho( 5 l7 -> 5 Is) mentioned earlier. The 
cross relaxation process is pictured in figure 7-4. 
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Figure 7-4 Illustration of the Cross-relaxation Process Ho( 5 l5-> 5 l7);Tm( 3 H6 -> 3 F 4 ) 
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The decay of the Ho 5 Is was measured for several co-doped samples with 
different Ho and Tm concentrations. Figure 7-5 shows a representative decay of the Ho 
5 l 5 at 912nm for excitation at 890nm. In this case Ho(.5)YLF and Tm(5)Ho(.5)YLF are 
pictured. As can be seen in this figure, the shortening of the lifetime is obvious. 

For each co-doped sample, the decay curve was fit to a series of exponentials 
using the Kaleidagraph™ data analysis/graphics application. It was found that a good fit 
was obtained with two exponentials. Table 7-1 tabulates the results of the fitting. The 
time-dependent fluorescence of the co-doped materials shows an initial quick exponential 
decay (~ 2 - 4 |xs) followed by a slightly longer exponential decay (-11-14 jis). 

We would like to be able to extract a rate constant for this process. Toward this 
end, we take a rate equation approach and write the rate of change of the Ho 5 Is 
population as. 


dn< nc 

“df = -^“P51 n 5 n l +P27 n 2«7 


(7-1) 


Table 7-1 Components of Ho 5 l 5 Decay @912 nm for Ho 5 Is Excitation @890 nm 


V. /> : y : : ;/••• 


n(t)/n(0) = Aiexp(-t/Ti) + A 2 exp(-t/T 2 ) 


Sample 

Ai 

Ti 

a 2 

'll 

Integrated 
time (fis) 


1.0000 

20.0 ! 

— 

— 

20.0 

EHS5HSJES 

0.5698 

2.67 

0.4543 

13.16 

7.47 

ESBBSEBB 

0.3938 

1.71 

0.6338 

11.92 

8.25 

HSSBBS 

0.4227 

2.75 

0.5788 

11.81 

7.96 

WSSBSSBS 

0.5587 

3.89 

0.4329 

13.43 

7.94 

1 SSSMiMS 

0.6294 

3.61 

0.3540 

11.58 

5.98 


0.4753 

66.5 

0.4712 

518 

280 
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Figure 7-5 Representative Decay of Ho 5 Is under Ho 5 Is excitation. 
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The subscript numbers refer to a particular manifold as indicated in figure 7-4. The 
first term in this equation refers to the natural decay of the Ho 5 Is manifold. The other 
terms account for the interaction of the manifolds participating in the process. In this case 
p51 is the rate constant of the cross relaxation Ho( 5 l5~> 5 l7);Tm( 3 H6 -> 3 F4) and p27 is the 
rate constant of the upconversion Tm( 3 F4-> 3 H6);Ho( 5 l7 -> 5 l5). Considering only the 
manifolds participating in these processes, we can write the ground state populations as 
nj = CtmN - n2, and ns = ChoN - ns - n7 . This is just a statement of conservation of Tm 
and Ho ions. Cxm is the concentration of Tm ions and N is the density of Y 3+ sites in 
YLF. Using these expressions in eq. 7-1 and neglecting all nonlinear terms, that is terms 
which are a product of populations, we find that: 

Y = -^-P51 c TmNn 5 (7-2) 

Neglecting the non-linear terms is a rather drastic assumption, but it is one which 
may not be too bad considering that the manifold which is being excited in the situation 
we are studying is the Ho 5 Is. It will initially have a large population relative to the other 
manifolds in the system and the non-linear terms will be unimportant. At later times, 
however, this assumption could drastically breakdown as the other manifolds also become 
populated. We will assume that the assumption holds fairly well and proceed. Equation 
7-2 is easily solved. The solution is: 

^ = Ae*p[-£-p 51 C Tnl N,] (7-3) 

This equation was fit to several Ho 5 Is decay curves with psi as the fitting 
parameter. T5 is known from measurements of singly doped YLF:Ho. The value of X5 was 
given in table 7-1 as 20ps. The results of the fitting appear in table 7-2 
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Table 7-2 Rate Constants For Cross Relaxation Ho( 5 l5-> 5 l7);Tm( 3 H6 -> 3 F4) 



n 5(t) _ 
n 5 (0) 


A exp 


'^■ _ P51 C Tm Nt 


HESS359MI 

A 



0.8598 

1. 4043xl0- 22 


0.8326 

1.0322x10-22 


0.8753 

1.2077x10-22 


0.8973 

1.0403x10-22 


0.8680 

1.3449x10-22 

Average 

0.8666 

1.2058x10-22 


It should be pointed out that the cross-relaxation process discussed here is not the 
only possible explanation of the Ho 5 Is decay in co-doped systems. It is only a possible 
explanation supported by the fact that the energy transfer is resonant, that is, the 
participating levels overlap in energy (see tables 5-1 and 6-1 for range of transitions ), and 
by the fact that the decay in Ho is exponential, eliminating an explanation based on 
interactions between Holmium ions in other states (see figure 5-7). The opposite process 
to the one described here, that is the upconversion Tm( 3 F4-> 3 H6);Ho( 5 l7 -> 5 Is), will be 
discussed next. 

7.2.2 Upconversion Tm( 3 F 4 -> 3 H 6 );Ho( 5 l7 -> 5 I 5 ) in YLF:Tm,Ho 

In section 7-1 we showed that there are several upconversion processes which 
may affect the laser performance of the 2pm laser in Tm,Ho systems. The upconversion 
process which may have the most deleterious affect on these systems has been identified 
as Tm( 3 F4-> 3 H6);Ho( 5 l7-> 5 l5), in which an excited Tm 3 F 4 ion interacts with an excited 
Ho 5 l7 ion. The Tm 3 F 4 ion decays back to the ground state and the Ho 5 l7 ion is excited 
to the Ho 5 l5. This process is illustrated in figure 7-6 The excitation in the 5 Is then decays 
back to the lower manifolds through the emission of phonons with a quantum efficiency 
that is near one. 
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Figure 7-6 Illustration of the Upconversion Process Tm( 3 F4-> 3 H6);Ho( 5 l7-> 5 l5) 


We would like to be able to measure the rate of this upconversion. Toward this 
end, a simplified rate equation model can be utilized. It will be useful to make some 
observations regarding the Tm,Ho system, when exciting the Tm 3 H4 manifold, to 
simplify the analytic equations governing the rate of change of the populations in the 
lower manifolds. One observation can be made from measurements of the upper state 
lifetimes in the Tm 3 F4 and Ho 5 l7 manifolds. Regardless of the Tm and Ho 
concentrations the upper state lifetimes in these two manifolds are the same. The 
lifetimes of these manifolds are very long compared to the time for the energy transfer 
from the Tm 3 F4 to the Ho 5 l7, so a state of quasi -thermal equilibrium is established 
between these manifolds. This results in a great simplification of the problem because the 
two upper states can be treated as a coupled system, and the relative populations of each 
Tm and Ho ion can be calculated from a Boltzmann distribution over the individual Stark 
levels. The fractional population fi m in the Tm 3 F4 and^j 0 in the Ho 5 l7 can be 
expressed as [Ca 75 ]: 
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(7-5) 


/Ho = 


~ v «- E j/ kT 

j 


C TmXSi e E * /kT + C HoSgj e Ej/kT 


The ion concentration of Ho is Cho and the ion concentration of Tm is Cx m The 
sums over i and j are over Stark levels of the Tm 3 F 4 and Ho s l 7 , respectively. These 
sums are the partition functions for the manifolds, and have already been calculated (see 
table 5-3 and table 6-3). The degeneracy of a given Stark level is taken into account by gi 
and gj. The energy of the Stark levels are Ej and Ej. 

Since the Ho 5 l 7 and the Tm 3 F 4 upper states can be treated as a coupled system, 
we can write down the rate equation for such a coupled system as [Fa88]: 


^ar = R P - - qN H oN T m + q \ N Ho N Tm cm) 


where N u =Nho +Nxm represents the population of the coupled system and X represents 
the coupled lifetime (the lifetime of the Ho 5 l 7 or the Tm 3 F 4 ). Now, this balance 
equation can be understood in the following way. The rate of change of the population in 
the coupled system is equal to the pump rate Rp (number of photons per second entering 
the coupled system within the pump volume of the crystal), minus the excitations lost 
through spontaneous emission (the 2nd term), minus the excitation lost through 
upconversion from the interaction of the Tm and Ho upper states (the 3rd term 
characterized by a an upconversion coefficient q), plus the excitations that return from 
decay of the 5 Is (the 4th term which includes T] a as the quantum efficiency for the 
relaxation to return to the Ho 5 l 7 and Tm 3 F 4 ). 

Now, we can write Nho and Nxm in terms of N u and the fractional populations 
fuo and /rm Since N u =Nho +Nxm then it follows that Nho^/HoNu and Nxm^/TmN u . The 
balance equation 7-6 now becomes very simple: 
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dN u = n _NjL_ n fi_!kl/ f N 2 
dt R P X 2 J-'Ho/TmNu 


(7-7) 


For a system in steady state this equation is easily solved. The result is: 

<7 - 8) 


A series of experiments have been performed on a number of YLF:Tm,Ho ciystals 
that measure the change in intensity of the Ho 5 l 7 intensity with absorbed pump power. 
A 250 mW diode laser operating at 785 nm was focused to a spot of 70 |xm x 130 Jim 
radius at the crystal face. The power of the diode laser was varied from 100%-5% of the 
total power with a variable attenuator. The absorbed power was measured by subtracting 
from the input power the power measured after passage of the diode laser through the 
crystal. Since the Ho 5 l 7 intensity is proportional to the population, then eq. 7-8 can be 
fit to the measured intensity vs. pump power curves to find the upconversion rate 
constant q. However, the population will depend on the spatial distribution of the pump 
volume inside the crystal. R p can be written in terms of the total pump rate R and an 
elliptical pump beam distribution [Fa88]: 


D _ ^P^ abs 

p " hv p 


2a 


jcw x w y [l - exp(-a£)] 


exp[-az] • exp 



cs 

M 



X , X 


z. 

2 + 2 






(7-9) 


In this equation Pabs represents the total absorbed power, hvp is the energy of the 
pump photons, Tip is the pump quantum efficiency, a is the absorption coefficient at the 
pump wavelength, w x and w y are the elliptical radii of the pump beam at the crystal face, 
and i is the thickness of the sample (path length). 
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Using eq. 7-9 in eq. 7-8 and integrating the population over the pump volume 
gives a value proportional to the intensity. All that needs to be done is to fit the intensity 
data to find the upconversion coefficient q. We did this for six samples of various Tm,Ho 
concentrations (five for YLF and one for YAG). The fits of eq. 7-8 to the data are shown 
in figures 7-7a-7-7f. The results are summarized in table 7-3. It is clear from this data that 
for moderate concentrations the upconversion coefficient is about an order of magnitude 
smaller in YLF than in YAG. 

We tried to determine what the concentration dependence of the upconversion rate 
might be. It was expected that the rate should depend linearly on the product of the 
Thulium and Holmium concentrations, but no such dependency could be determined from 
the experimental results. The q values appear to follow a linear dependence with the ratio 
of the concentrations, but this makes little physical sense. 


Table 7-3 Parameters Used in the Data Fit to Upconversion Coefficient q 


Sample 

x(ms) 

fia 

/ho 

/Tm 



16.5 ms 

1.91 

0.2940 

0.7060 

7.0x10*24 


19.0 ms 

1.91 

0.5554 

0.4446 

2.3x10*24 


15.0 ms 

1.95 

0.1176 

0.8824 

1.0x10*23 


13.0 ms 

1.97 

0.3570 

0.6430 

6.5x10-24 


8.5 ms 

2.00 

0.5260 

0.9473 

4.0x10-23 

EBSSS&S59 

12.0 ms 

1.95 

0.4179 

0.5821 

1.2x10-23 
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Figure 7-7a Fitting for Upconversion Coefficient in Tm(4)Ho(.5) YLF 
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Figure 7-7b Fitting for Upconversion Coefficient in Tm(4)Ho( 1 . 5)YLF 
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Figure 7- 7c Fitting for Upconversion Coefficient in Tm(5)Ho(.2) YLF 
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Figure 7-7d Fitting for Upconversion Coefficient in Tm(12)Ho(.2) YLF 
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Figure 7-7e Fitting for Upconversion Coefficient in Tm(6)Ho(l) YLF 
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Figure 7-7f Fitting for Upconversion Coefficient in Tm(6)Ho(l) YAG 
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It should be noted that q is identified with the rate constant P27 in the notation 
used in section 7.2.1, representing the opposite process described by ps|. As can be seen 
from the values tabulated in table 7-3, there is a bit of scatter in the values for the samples 
studied. If we take an average of the YLF values, neglecting the sample with 12% Tm we 
find that q = P27 = 6.45x1 0* 24 cm 3 4is. This value should be considered, at best, a rough 
estimate of the rate constant due to the convoluted way in which it was determined. One 
of the major problems, aside from the assumptions of the rate equation model, is the 
measurement of the excitation volume. If the measured volume is in error by a factor of 2, 
the q value could change by an order of magnitude. A better way to determine p27 would 
be to pump the Tm 3 F 4 manifold directly and look at the decay from the Ho 5 l5. 
Unfortunately, no signal could be detected for this emission under our experimental 
conditions. 


7.2.3 Cross-relaxation Tra( 3 H 4 -> 3 H6);Ho( 5 I 7 -> s S 2 ) in YLF:Tm,Ho 

The decay of the Ho 5 S2 exhibits a similar behavior to the Ho 5 Is, that is, the 
decay is exponential for singly doped Holmium, but shortens and becomes non- 
exponential with the addition of Thulium. Figure 7-8 shows the decay curves for the 5 S 2 
manifold when exciting the s S2 manifold in co-doped systems with a dye laser. The data 
for these decays were received courtesy of G. Armagan. The behavior is non exponential 
with an initial quick component (~ 9|is) followed by a longer one (~ 55|is). This is 
interpreted as the cross-relaxation process Ho( 5 S2-> 5 l7);Tm( 3 H6 -> 3 H4). We can make a 
similar rate equation analysis as was done for the 5 Is cross-relaxation. In this situation it 
was found that the rate constant for the cross-relaxation Ho( 5 S2-> 5 l7);Tm( 3 H6 -> 3 H 4 ) 
was p9i = 2.59xl0' 23 cm 3 /ps. In the notation here, manifold 9 corresponds to the Ho 5 S2 
and manifold 1 corresponds to the Tm 3 H6- 
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Figure 7-8 Luminescence Decay of 5 S 2 Manifold After Excitation of 5 S 2 @535 nm 
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7.2.4 Energy Transfer Tm( 3 F 4 -> 3 H 6 );Ho( 5 l8-> 5 l7) in YLF:Tm,Ho 


We now turn to a discussion of the energy transfer between the Tm 3 F 4 and Ho 
5 l 7 manifolds in co-doped YLF. The transfer of energy from the Tm 3 F 4 to the Ho is 
of primary interest since it is by this process that the Ho 5 l 7 manifold becomes 
populated, and makes possible laser action at 2\m. The transfer has been interpreted by 
the process Tm( 3 F 4 -> 3 H 6 );Ho(%-> 5 l 7 ) and is pictured in figure 7-9. 


3 H 4 

3 h 5 

3 F< 

3 h 6 
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® 
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© — 

© T 
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5 1 5 
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5 1 7 
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Tm 3+ Ho 3+ 

Figure 7-9 Illustration of the Energy Transfer Process Tm( 3 F 4 -> 3 H 6 );Ho( 5 l 8 -> 5 l 7 ) 


This process is a resonant one since the participating levels exhibit very good 
overlap (see tables 5-1 and 6-1) and is supported by the decay characteristics of these 
two manifolds. Figure 7-10a and 7-10b shows the time dependent fluorescence of the Tm 
3 F 4 and Ho 5 l 7 respectively, for two different excitation wavelengths (a TiAl 2 C >3 tuned to 
785nm and a CoMgF 2 tuned to 1.7|xm). Figures 7-10 illustrate nicely the correspondence 
between the decay of the Tm 3 F 4 and the rise of the Ho 5 l 7 . For the purpose of studying 
this energy transfer process, direct excitation of the Tm 3 F 4 greatly simplifies the analysis 
of the decay curves. For example, in figure 7- 10a, when exciting the 3 H 4 , the 3 F 4 decay 
has a rise time which reflects the buildup of the population as a result of the cross- 
relaxation Tm( 3 H 4 -> 3 F 4 );Tm( 3 H6-> 3 F 4 ) and the natural decay of the 3 H 4 . 
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Normalized Intensity 



Figure 7-10a Decay of Tm 3 F 4 in Tm(4)Ho(.5)YLF For Different Excitation Wavelengths 



Figure 7-10b Decay of Ho 5 l 7 in Tm(4)Ho(.5)YLF for Different Excitation Wavelengths 
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This complicates the analysis because the mixture of many processes combine in 
the decay profile and it becomes difficult to distinguish between various components. On 
the other hand, by exciting the Tm 3 F 4 a much simpler decay profile develops without the 
complication of the rise time. However, since most of the research on Tm,Ho lasers have 
utilized Alexandrite or diode pumping of the Tm 3 H 4 to achieve laser action, the time 
dependent fluorescence for Tm 3 H 4 excitation is perhaps more useful. From the 
spectroscopic point of view, however, we will concentrate mainly on the decay 
characteristics of the Tm and Ho manifolds for Tm 3 F 4 excitation. 

In figure 7- 10b we see an obvious difference in the buildup of the population in 
the Ho 5 l 7 when exciting the Tm 3 F 4 and the Tm 3 H 4 . The rate at which the Ho 5 l 7 is 
populated is slower for Tm 3 H 4 pumping than for Tm 3 F 4 pumping. The reason for this 
is that for Tm 3 H 4 pumping the Tm 3 F 4 is populated mainly via cross-relaxation and 
increases gradually with time. During this time energy transfer is also occurring between 
Tm and Ho, but is dependent on the number of excited Tm ions available for transfer. So, 
the buildup in the Ho 5 l 7 is dependent on the efficiency and rate of the cross-relaxation 
process. On the other hand, directly pumping the Tm 3 F 4 creates a large population of 
Tm ions available for transfer. As a result, the Ho population rises much faster. 

We will now examine the details of the Tm -> Ho energy transfer process 
Tm( 3 F 4 -> 3 H 6 );Ho( 5 l 8 -> 5 l 7 ). We have said before that this is a resonant process since the 
energy levels of the Tm 3 F 4 , 3 H6 and Ho 5 l 7 , 5 Is overlap. This is shown pictorially in an 
energy level diagram in figure 7-11. The overlap is quite good and considering that the line 
width for transitions at room temperature is about 15-20 cm' 1 there are a large number of 
resonant transitions possible between Stark split levels in the lower Tm and Ho 
manifolds. For an energy mismatch of only 1 cm' 1 there are three transitions possible: 

1. ) Tm(5599 ->305);Ho(0 ->5293) 

2. )Tm(5599 ->372);Ho(7 ->5233) 

3. ) Tm(5599 ->419);Ho(48 ->5229) 
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In general, the process can be written as Tm(i ->j);Ho(f ->k) where i, j, £, k 
denotes a Stark level of the 3 F 4 , 3 H6 , 5 Is , and 5 l 7 manifold respectively. The probability 
of finding a Tm 3 F 4 ion in Stark level i is given by the Boltzmann distribution Pj 


exp(-AEj/kT) 
Pi £exp(-AEi/kT) 


(7-10a) 


and the probability of finding a Ho 5 Ig ion in Stark level l is given by: 

exp(-AEi / kT) 

Vl £exp(-AE*/kT) 

t 


(7- 10b) 


Note that AEj or AE^ are energy differences from the lowest Stark level to the Stark 

level i or l in the 3 F 4 or 5 Is manifold. So, the energy transfer rate for a given transition 
Tm(i ->j);Ho( £ ->k) will depend on the product of the probability of an ion being in Stark 
level i of the Tm 3 F 4 and an ion being in Stark level £ of the Ho 5 Is and a temperature 
independent transfer rate constant, which will contain matrix elements between the initial 
and final states. This can be written as: 


hi = Wtf p^ 


(7-1 la) 


Actually, there are many possible transfers both resonant and non resonant, so eqs. 
7-10 must be summed over all transitions that are deemed physically realizable. The rate of 
backtransfer will have a similar form: 


^kj = w kjPkPj 


(7-llb) 
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Figure 7-11 Stark Levels of the Ground and Lower Manifolds in YLF:Tm,Ho 
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It is obvious that such an approach to calculating energy transferrates is intractable 
given the number of possible transitions, assuming one can even identify all of them. 
Nevertheless, the exercise is useful in understanding how energy is transferred. It also 
shows us that describing the population of the combined Tm 3 F4 and Ho populations 
by a thermal distribution is justified. That this is true follows from the large number of 
possible transitions and the fact that the transfer times are much faster than the lifetimes of 
the manifolds. At any rate, we can write down the system of equations of the rate of change 
of the Tm 3 F4 and Ho 5 l7 populations. 

The rate equations are: 


dn2 _ n2 
dt t 2 


- P28 n 2 n 8 + P71 n 7 n l 


( 7 - 12 a) 


dni n7 

_ = + P28 n 2 n 8 - p 71 n7n! 


( 7 - 12 b) 


In these equations the subscript numbers refer to a particular manifold as denoted in 
figure 7 - 9 . %2 and T7 are the natural (radiative + nonradiative) lifetimes of the Tm 3 F4 and 
Ho 5 l7 manifolds the absence of energy transfer. P28 and P71 are the rate constants of Tm 
-> Ho transfer and Ho -> Tm backtransfer respectively. 

Now, if we assume that these are the only manifolds with any sufficient 
population, which is probably the case for direct excitation of the Tm 3 F4 , then we can 
write ni + n2 = CjmN and n7 + ns = ChoN, where Cx m and Cho are the dopant 
concentrations of Tm and Ho. N is the number of Y 3+ sites in YLF into which the dopant 
ion are substituted. This allows us to eliminate the ground state populations from the rate 
equations and write them as: 


(in 2 _ n2 
dt ^2 


- otn 2 + pn 7 + (p2g - P7i)n 2 n 7 


( 7 - 13 a) 
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(7-13b) 


__ .« + an 2 - Pn 7 - (p 2 8 - p 7i )n 2 n 7 

where we have written a = P 28 QhoN and p — p 7 iCTmN. a and P are the rates of energy 
transfer and back transfer analogous to eq. 7-1 la and 7-1 lb. This is the reason why we 
referred to P 2 S and p 7 i as rate constants because they are independent of concentration. 
The set of coupled equations can be solved numerically for supplied values of X 2 . T 7 , a 
and p. We have written a computer program which solves these differential equations 
numerically that uses a Runge Kutta method. The solutions are sufficiently smooth that it 
was not necessary to incorporate a variable step size for differentiation at various points. 
It was discovered that the solutions did not vary much if the nonlinear term involving 
n 2 n 7 was dropped, as long as the initial condition of the ions excited into the upper 
manifold at t=0 were relatively small. This condition corresponds to low pump densities 
of the excitation source. For the CoMgF 2 laser we used only had 6mj of energy in a beam 
size that was 100’s of microns in diameter. This allows us to linearize equations 7-13 and 
solve them to obtain closed form solutions for the rate of change of the populations n 2 and 
n 7 . The solutions were found to be (the derivation is given in Appendix D): 

^=(^) exp( - ,/t)+ (^f) xp[ " (a+p)t) (7 ' 14a) 

^ + P >1] < 7 - 14b) 

The coupled lifetime x in equations 7-14 represents the lifetime of the Tm and Ho 
manifolds since after a short amount of time the energy levels in Tm and Ho reach a thermal 
equilibrium and decay at the same rate. This is indeed what is found from the experiments. 
After several hundreds of microseconds, the decay of Tm and Ho proceed at the same rate. 
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This is due to the nature of the Tm <-> Ho sharing process. The upper manifolds have a 
certain amount of energy as a combined system and an agreement through Boltzmann 
statistics to maintain the distribution between the two. Even though the total amount of 
energy in the upper states is decreasing, the distribution between Tm and Ho is maintained, 
so they decay at the same rate. By looking at the very early parts of the decay curves, 
however, we can catch Tm and Ho in the act of sharing their energy. This is exactly the 
energy transfer process we are interested in. 

Decay Measurements using the CoMgF2 laser tuned to 1.75jim the were made on 
several co-doped YLF samples containing various concentrations of Tm and Ho. The 
emission was monitored at 1.7pm and 2.06pm and the decay curves obtained. The 
intensity was normalized to one and the decay curves for the Tm 3 F4 were fit to equation 7- 
14a to obtain the energy transfer rates a and p. The values of a, P and T obtained by 
fitting the decay curves to Eq. 7-14a are shown in table 7-4. Figures 7-12a - 7-12h 
display the first 1000 ps of the measured decay curves of eight YLF samples co-doped 
with different amounts of Tm and Ho. The best fit to the data is also displayed in these 
figures by the dashed lines. 

Now that we have the rates a and P , we can extract the rate constants and see 
how well they agree with each other because, as we said earlier, they should be 
independent of concentration. Recall that we defined a = P28QloN and P = P7iCTmN. 
We know the dopant concentrations and the density of sites in the samples. For YLF and 
YAG N=1.38xl0' 22 ions/cm 3 and for LuAG N=1.42xl0' 22 . It is therefor a simple matter 
to extract the rate constants P28 and P71 . The results are presented in table 7-5. It is seen 
from these calculated rate constants that they do not vary too much for the different 
samples. Overall, the agreement is not too bad considering that the fits to the data were 
marginally okay. 
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Figure 7-12h Decay of Tm(l)Ho(12)YLF 
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Table 7-4 Values of the Fitting Parameters to eq. 7-14a 


■Knxi 

a (Rs‘ l ) 

P (US' 1 ) 

X(ns) 


0.004125 

0.004532 

15858 

BSESSfflSBSI 

0.009176 

0.005243 

16237 


0.012221 

0.007178 

18674 

ESS5HEB33 

0.002750 

0.006295 

14502 


0.005863 

0.007103 

10491 


0.009152 

0.011415 

13103 

B93S393S39 

0.005922 

0.010284 

25108 


0.020101 

0.002016 

8519 


0.005124 

0.011145 

11758 


0.002901 

0.005944 

9892 


Table 7-5 Values of the Rate Constants for Energy Transfer Tm( 3 F 4 -> 3 H 6 );Ho( 5 l 8 -> 5 l 7 ) 


Sample 


K9E15HI 

WMSBSHM 

IWWBBH 

0.5978x10-22 

0.0821x10*22 

0.1373 

BHBWBB31 

0 6649x10-22 

0.0949x10-22 

0.1427 


0.5898x10-22 

0.1300x10*22 

0.2204 


0.9963x10-22 

0.0912x10-22 

0.0915 

Tm(5)Ho(0.5)YLF 

0.8497x10-22 

0.1209x10-22 

0.1422 

EB5SS9SBEI 

0 6631x10-22 

0.1378x10-22 

0.2078 

BHBSgffiSSHI 

0.4913x10-22 

0.0621x10-22 

0.1263 

BHBSSSBBH 

0.7304x10-22 

0.1268x10-22 

0.1736 

(H55BES0 

0.3713x10-22 

0.1346x10*22 

0.3625 

Tm(5)Ho(.5)LuAG 

0.4083x10-22 

0.0836x10-22 

0.2047 
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The Tm->Ho transfer and Ho->Tm transfer interaction parameters can be found 
by calculating the appropriate spectral overlap integral in the Forster-Dexter theory of 
resonant energy transfer. Using eq. 6-1 we found that: 

Tm->Ho transfer Tm( 3 F 4 -> 3 H 6 );Ho( 5 l8-> 5 l7): C D a= 1.6901X10- 45 cm 6 /ns 
Ho->Tm transfer Ho( 5 l7-> 5 l8);Tm( 3 H6-> 3 F4): C D a= 1.2424x1 O' 46 cm 6 /ps 
Tm->Tm migration Tm( 3 F4-> 3 H6);Tm( 3 H6-> 3 F4): Cdd = 2.6191x10' 45 cm 6 /tis 
Ho->Ho migration Ho( 5 l7-> 5 l8);Ho( 5 l8-> 5 l7). Cdd = 7.2107x1(H 5 cm 6 /|i.s 

The connection between the microscopic interaction parameters and the rate 
constants has not been put into a theoretical framework at this time, however, we would 
expect that the ratio of the microscopic interaction parameters should be the same as the 
ratio of the rate constants for Tm->Ho transfer and Ho->Tm transfer. The ratio of the 
interaction parameters in this case is CDA(Ho->Tm)/CDA(Tm->Ho) = 0.0735. This is a 
factor of two smaller than the average value of the ratio of the rate constants in table 7-5, 
where it is found that p7i/p28 = 0. 1552. 

The reason for the difference in these two ratios possibly resides in the fact that 
the ratio calculated with the interaction parameters accounts only for the resonant energy 
transfer via the spectral overlap, whereas the ratio of the rate constants accounts for all 
transfers both resonant and non-resonant. To see if this reasoning is correct, we have 
derived a relation in terms of known quantities for the ratio of the rate constants assuming 
only the resonant transitions participate in the process. The derivation is given in 
Appendix E. The result is: 

5£=ifH-( E * 4 - E ™)/ k B T l <7 - l6 > 


131 



where Z\ is the partition function of the Tm 3 H 6 , Z 2 is the partition function of 
the Tm 3 F 4 , Z 7 is the partition function of the Ho 5 l 7 , and Zg is the partition function of 
the Ho 5 l 8 . and E ^ 7 are the energies of the lowest Stark levels in the Tm 3 F 4 and 
Ho 5 l 7 manifolds, respectively. The values for all of these terms can be found in table 5-3 
and table 6-3. The ratio of the rate constants for T=300“K was calculated from eq. 7-16 
and found to be p7i/p28 = 0.0616. This is very close to the ratio of the interaction 
parameters, giving credence to our earlier statement. 

We can use the ratio of Tm->Ho transfer and Ho-> Tm backtransfer to 
experimentally determine the fraction of ions residing in the excited Ho ions after 
equilibrium between the Tm and Ho ions is reached. That is: 

/Ho =7^ < 7 ' 17 > 


which can be rewritten as: 


/Ho = 


1 

l + (n 2 /n 7 ) 


(7-18) 


In steady state, we can set dn 2 /dt = dn?/dt = 0. That is we can set equations 7-12 
equal to zero. If we set eq. 7- 12a equal to zero and neglect the term n 2 /t 2 as being small 
compared to the others, then we obtain. 


P 71 _ n 7 n I 
P28 n 2 n 8 


(7-19) 


Using ni = n 2 - Or m N and ns = n 7 - ChoN, the expressions for conservation of 
Tm and Ho ions, and neglecting the non-linear terms involving n 2 n 7 we find that: 
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This can be used in eq. 7-18 to find the fraction of Ho ions in the state. The 
results are tabulated in table 7-6 

Table 7-6 Fraction of Ho Ions in the 5 l 7 state For a Given Concentration 



n2/ii7 

Jib) 

II 

i—i 

t 

BSBSSJSSBBil 

1.2416 

0.4461 

0.5539 

WSGBM5SM 

0.6208 

0.6169 

0.3831 


0.4138 

0.7072 

0.2928 


3.8800 

0.2049 

0.7951 

EfSSSH 

1.5520 

0.3918 

0.6082 


0.9312 

0.5178 

0.4822 

E5555BS3I 

0.0129 

0.0969 

0.9031 


9.3120 

0.9872 

0.0128 


Of course, this is not the only interpretation of the data that can be put forth. In fact, 
given that the rate equation method did not fit the data better would suggest that the decay is 
not composed of two exponentials, but perhaps three or more. We did a curve fit using 
Kaleidagraph™ data analysis/graphics application to fit the data to three exponentials. 
Kaleidagraph™ uses the Levenberg-Marquart algorithm to fit the data. 

Six parameters were used that composed an equation for the sum of exponentials of 
the fonn: 


= A! exp(-t /t 1 ) + A 2 exp(-t / x 2 ) + A 3 exp(-t / x 3 ) 


(7-21) 


























The best set of parameters Aj, ti, A2, t2,A3, 13 were obtained from the curve fits 
of the experimental data performed by Kaleidagraph™. The data fits were excellent for 
three exponentials. The parameters appear in tables 7-7 and 7-8. 

From the data in tables 7-7 and 7-8 we can see first of all that the long time decay 
of Tm and Ho are the same. This is as expected because of the thermalization between 
Tm and Ho discussed earlier. The early time behavior is quite different however. In both 
cases the early time behavior (0 - 1000 ps) appears to be composed of two exponentials 
with time constants in the tens to hundreds of microseconds. The rise time of the Ho 5 l7, 
however, seems to be faster than the initial decay time of the Tm 3 H4- This would be 
contrary to common sense because if Tm is feeding Ho initially before thermal 
equilibrium, then the decay in Tm would be expected to match the rise in Ho. It must be 
remembered that the signal detected that forms the decay curve is the fluorescence, which 
is due to radiative emission from the ions. Given a sufficient amount of time, an ion in Ho 
or Tm will emit radiatively unless it finds somewhere to transfer its energy or can give 
phonons to the lattice and de-excite that way. 

At t=0 Tm has many excitations and Ho has none. Immediately Tm starts 
transferring to Ho because all Ho sites welcome the energy. Not many of the Tm ions can 
emit radiatively because Ho is stealing the excitation energy from Tm. Many Ho ions 
begin to exist in excited states, but still less than Tm. The Ho ions distribute themselves 
among the Stark levels, but in levels which are less resonant with Tm. Remember that at 
the beginning of this section we saw that the most resonant transitions come from the 
bottom of the 3 F4 (5599 cm -1 ) and they end up at the top of the 5 l7 (5293 cm -1 , 5233 
cm* 1 , and 5229 cm -1 ). When the excitation distributes itself in the lower Stark levels of the 
Ho 5 l7 then backtransfer is less resonant. The reader can convince oneself by examining 
figure 7-11. So, while the backtransfer exists, it should be slower than the transfer. Thus, 
initially Ho keeps its excitation longer and increases its chances of decaying radiatively. 
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Table 7-7 Exponential Components of Tm 3 F 4 Decay for Tm 3 F 4 Excitation 


Sample j 

— 

Al 

RI81 

mm 

in 

/ .o. 

A 3 

a aa— 

dffil 

mSSSBESSM 

0.3411 

57.5 

0.1802 

319.8 

0.4721 

15998 


0.4656 

37.6 

0.1667 

265.2 

0.3591 

16474 


0.4451 

26.8 

0.1841 

165.3 

0.3677 

18412 

Tm(5)Ho(0.2)YLF 

0.1887 

66.2 

0.1138 

366.5 

0.6911 

14596 


0.3166 

47.7 

0.1359 

312.5 

0.5402 

10633 

W5SSSSEBSM 

0.2688 

23.1 

0.1735 

146.2 

0.5509 

13197 

Tm(l)Ho(12)YLF 

0.3960 

38.8 

0.2116 

195.2 

0.3625 

25310 


0.0835 

50.9 

— 

— 

0.9068 

8536 

Tm(6)Ho(1.0)YAG 

0.2256 

39.4 

0.0770 

354.0 

0.6762 

11907 

Tm(5)Ho(.5)LuAG 

0.1647 

48.9 

0.1622 

282.5 

0.6650 

9926 


Table 7-8 Exponential Components of Ho Decay for Tm 3 F 4 Excitation 


Sample 

Ai 


a 2 

ISRI 

A 3 


Tm(4)Ho(0.5)YLF 

0.4754 1 

21.1 

0.4306 

183.5 

1.0669 

17125 

Tm(4)Ho(1.0)YLF 

0.4325 

5.65 

0.5004 

61.5 

1.0396 

16361 


0.5291 

4.76 

0.4585 

39.3 

1.0264 

19707 

Tm(5)Ho(0.2)YLF 

0.4333 

22.5 

0.3723 

213.8 

1.0814 

15116 

Tm(5)Ho(0.5)YLF 

0.4415 

9.07 

0.5531 

79.4 

1.0645 

10640 


0.6109 

5.36 

0.3723 

44.5 

1.0292 

13212 


1.0613 

2.52 

— 

— 

1.0009 

25123 

Tm(l 2)Ho(0.2)YLF 

0.7491 

4.32 

0.2893 

30.3 

1.0238 

8597 

WSBSSSEBSi 

0.7201 

4.97 

0.2897 

57.9 

1.0325 

11535 

Tm(5)Ho(0.5)LuAG 

0.4049 

22.8 

0.3205 

183.5 

1.0766 

10222 
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The result of all this is that the rise in Ho should be faster than the decay of Tm. 
As time goes on things begin to equal out and thermal equilibrium is reached between the 
Tm and Ho ions. However, since the backtransfer from Ho to Tm is less effective than 
the transfer of Tm to Ho, the thermalized lifetime of the combined Tm,Ho system decays 
at the lifetime that the Ho ions are decaying at. This will not always be the case because 
of radiative trapping or impurities that quench the lifetime. 

Finally, we suggest another theory which does not necessarily negate the rate 
equation model but suggests a way to modify it. The fact that there are two exponentials 
in the early part of the decay suggests that there may be two distinct regions of decay. In 
one region a Tm ion sees and enhanced region of Ho ions and another region where Tm 
sees a depleted region of Ho ions. The rate equation model as was proposed here does not 
account for changes in the donor and acceptor environments. This offers an explanation 
for the excellent fit of the decay by three exponentials, but does not suggest a feasible 
way to fix the rate equation model. A fusion between these two models and a better 
characterization of the energy transfer process on the lines of that discussed at the 
beginning of this section would perhaps be a good starting point for further investigation. 
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CHAPTER 8 

Summary and Conclusions 


The results of the preceding investigation lead to a very interesting conclusion, 
which is here to be deduced 

Albert Einstein 

Does the Inertia of a body depend 

upon its energy content? 

Annalen der Physik, 17, 1905 

8.1 Summary of the Results 

A detailed analysis of the spectroscopy and excitation dynamics of Tm 3+ and 
Ho 3+ in Lithium Yttrium Fluoride has been done. Quantitative as well as qualitative 
approaches have been taken towards the interpretation of the experimental data. This is 
particularly evident in chapter 7, where the various energy transfer processes are 
discussed and interpreted. A purely quantitative approach is limiting in scope and does 
not always lead to new ideas. It can, in fact, be misleading especially where modeling is 
concerned. The construction of a model that fits the data well does not imply the model is 
correct. The physical reasoning behind it must be sound also. So, attempts have been 
made to discuss at length the reasoning behind the ideas presented. 

The introductory chapters (chapters 1-4) laid the groundwork for this study. A 
thorough introduction was presented in chapter 1 to familiarize the reader with the 
subject under study, and to give a sense of the direction this thesis would take. A brief 
preview of the contents and results were also presented to peak the interest of potential 
readers. The theoretical background was given in chapter 2. This is material which is well 
known. The classic textbooks of Di Bartolo, Wyboume, Judd and Dieke really have no 
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substitute, all written during the 1960's when laser spectroscopy was a new field. 
Chapter 3 covered the details of the YLF host lattice. While the properties of symmetry 
and group theoretical treatment of the ciystal field effects on ions in solid state materials 
is a fascinating subject in itself, only the details and topics necessary to an understanding 
of the phenomena discussed here were presented. Chapter 4 carefully detailed the 
experimental apparatus and techniques used in the data collection. None of the 
experiments were novel or new, but are well worn favorites. 

Chapters S through 7 constitute the presentation of the data and the analysis. 
They form, in a sense, a trilogy: The results and interpretation of YLF.Tm, The results 
and interpretation of YLF:Ho, and the results and interpretation of YLF:Tm,Ho. The 
outcome of a Judd-Ofelt analysis on YLF:Tm and YLF:Ho were given. Manifold to 
manifold radiative lifetiimes and branching ratios were determined. Comparison of 
emission and absorption cross sections provided an alternate method for extracting 
radiative lifetimes.The correspondence between the radiative lifetimes calculated by the 
two methods was satisfactory. They both gave results in the same order of magnitude 
and, in fact, differed by less than a factor of two. 

A limited presentation and discussion of the excitation dynamics of singly doped 
YLF:Ho was given at the end of chapter 5. It was very interesting to discover that the 
green luminescence from the Ho 5 S 2 was observed when exciting the Ho 5 Is and also when 
exciting the Ho 5 F 5 , both of which lie lower in energy than the Ho 5 S 2 . Some kind of 
upconversion in Ho is evident from this observation, but all of the manifolds studied had 
nearly the same decay behavior whether they were excited directly or not. The conclusion 
is that this upconversion is rather weak or the result of excited state absorption by a 
second pump photon. 

In chapter 6, the excitation dynamics of YLF:Tm was discussed. Under 3 H 4 
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excitation, the concentration dependence of the 3 H4 fluorescence was examined. The 
concentration quenching is generally interpreted as a process in which an ion initially 
excited to the 3 H4 state interacts with an unexcited ion in the 3 H$ ground state. The ion in 
the 3 H4 state gives up part of its energy to the ion in the 3 H6 state, producing two ions in 
the 3 F4 state.This energy transfer process is called cross relaxation and is referred to 
notationally as Tm( 3 H4-> 3 F4);Tm( 3 H6-> 3 F4). This process was seen to be accompanied 
by energy migration amonq Tm ions within the 3 H4 state Tm( 3 H4-> 3 H6);Tm( 3 H6-> 
3 H4). The migration was interpreted as a hopping process in which the energy migrates in 
a random walk process among ions in the 3 H4 state until it happens to hop onto a 
quenching center where cross relaxation takes place. The hopping model was used to find 
the interaction strengths of the migration and cross relaxation processes (Cda, and Cdd) 
by fitting the model equation to measured decay curves of the 3 H4 fluorescence. The 
values of die interaction parameters obtained from the hopping model were compared to 
the values obtained from the Forster-Dexter theory by calculating the appropriate 
spectral overlap integral. Good agreement was found between these two methods. In both 
cases it was found that Cda < Cdd» within the range of applicability of the hopping 
mechanism. 

In chapter 7, various energy transfer processes in co-doped YLF:Tm,Ho were 
examined. Under cw diode laser excitation of the Tm 3 H4, co-doped YLF:Tm,Ho was 
seen to exhibit fluoresence at many different wavelengths. The transitions corresponding 
to these wavelengths were seen to originate from both Tm and Ho ions, indicating the 
presence of energy transfer from Tm to Ho. Some of these processes were examined by 
directly exciting the manifolds involved in the energy transfer with a pulsed laser source 
tuned to the appropriate wavelength. 

Of particular interest, due to its importance for understanding 2.0pm YLF:Tm,Ho 
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lasers, is the energy transfer from the Tm 3 F4 to the Ho 5 l7 Tm( 3 F4-> 3 H6);Ho( 5 l8-> 
5 l7). This transfer was found to occur in both directions, with the Tm->Ho transfer 
occuring more rapidly than the Ho->Tm transfer. A state of equillibrium between Tm 
ions in the 3 F4 and Ho ions in the 5 l7 sets in quickly ( 100 - 400 ps) after which the Tm and 
Ho ions decay at the same rate 

To probe into the dynamics of this energy transfer process the Tm 3 F4 manifold 
was excited directly with a CoMgF2 pulsed laser, and die time dependence of the Tm 3 F4 
and Ho 5 l7 fluorescence decay was measured. A rate equation model for this Tm->Ho 
energy transfer was developed, and used to model the decay behavior. For the excitation 
density in our experiments, it was discovered that the non-linear terms in the rate 
equations were unimportant. This allowed us to solve the set of coupled differential 
equations in closed form. The closed form solutions were fit to the experimental data 
collected for various Tm,Ho dopant ion concentrations. Rate constants for the Tm->Ho 
energy transfer and Ho->Tm backtransfer were extracted from the best fit to the data. The 
rate constants showed only minor variation within the range of Tm,Ho concentrations 
studied, indicating the success of the model, despite the fact that the model did not 
exactly fit the experimental data. The basis for this reasoning is that the rate constants 
represent an average transfer rate, and are only capable of predicting an average bahavior 
within the framework of the model. In reality, the transfer rates will depend on the spatial 
distribution of Tm and Ho ions in the crystal, and on the thermal distribution of 
excitations within the Stark levels of the manifolds involved in the energy transfer. 

Tables 8-1 and 8-2 compile some values generated by this thesis that may be 
useful for researchers interested in the design of YLF:Tm,Ho lasers. 
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Table 8-1 Useful Values and Parameters for YLF Tm,Ho 


Symbol and value 

Judd Qfelt Parameters Ho YLF 

Q 2 = 0.7483x10-20 cm 2 

0.2 - 1.6168x 10‘20 cm2 

0.2 = 1.2269x 10"20 cm2 

Judd Qfelt Parameters Tm YLF 

Q 2 = 1-6673x10-20 cm2 
Q 2 = 1.9184x10-20 cm2 
Q 2 = 1.0897x10-20 cm 2 

Energy level range in Ho YLF 

AE 517 (5153-5293cm-l), AE 518 (0-315cm-!) 

Energy level range in Tm YLF 

AE 3F 4 (5599-5972cm - l), AE 3H 6 (0-419cm*l) 

Partition functions in Ho YLF 

Z 517 = 11.415, Z 518 = 9.199 

Partition functions in Tm YLF 

Z 3 p 4 = 3.432, Z 3 h 6 = 5.207 

Radiative lifetimes in Ho YLF 

T 515 = 7.0 ms, T 516 = 5.6 ms, T 517 = 14.0ms 

Radiative lifetimes in Tm YLF 

T 3 h 4 = 2.4 ms, T 3 f 4 = 15.6ms 

Emission cross sections of Ho 5 l 7 

a = 1 .2286x10-20 cm 2 (n-pol.,E^c) @2.06^im 
a « 0.6530x10-20 cm 2 (o-pol.,E 1 c)@2.06*im 

Absorption cross section of 
Tm 3 H 4 @785 nm 

0.312x10-20 cm 2 (same for it-pol. and a-pol.) 





















Table 8-2 Measured Energy Transfer Parameters 


Value 

Interaction parameter of Tm-Tm cross- 
relaxation: Tm( 3 H 4 -> 3 F 4 );Tm( 3 H 6 -> 3 F 4 ) 

C DA - 4.7402x1 0" 47 cm 6 /ps 

Interaction parameter of Tm-Tm migration 
Tm( 3 H4-> 3 H6);Tm( 3 H6-> 3 H4) 

Cdd = 2.7238X10- 45 cmfyps 

Interaction parameter of Tm->Ho energy 
transfer: Tm( 3 F 4 -> 3 H 6 );Ho( 5 l 8 -> 5 l 7 > 

Cda = 1.6905x1 O' 45 cm 6 /ps 

Interaction parameter of Ho->Tm energy 
transfer: Ho( 5 l 7 -> 5 Ig)Tm( 3 H 6 -> 3 F 4 ) 

C D a = 1.2424xl0‘ 46 cm 6 /p.s 

Rate constant of Tm->Ho energy transfer: 
Tm( 3 F 4 -> 3 H6);Ho(5l 8 ->5i 7 ) 

P28 = 6.979xl0 -23 cm 3 /p.s 

Rate constant of Ho->Tm energy transfer: 
Ho( 5 l 7 -> 5 l 8 );Tm( 3 H 6 -> 3 F 4 ); 

P71 = 1.057x1 O' 23 cm 3 /|is 

Rate constant of Tm-Ho upconversion 
Tm( 3 F 4 -> 3 H6);Ho(5l 7 ->5i 5 ) 

P27 = 7.560xl0' 24 cm 3 /ps 

Rate constant of Ho-Tm cross relaxation: 
Ho( 5 I 5 -> 5 l7)Tm( 3 H6-> 3 F 4 ); 

P51 = 1.2058xl0' 22 cm 3 /|is 

Rate constant of Tm-Ho upconversion 
Tm( 3 H 4 -> 3 H 6 );Ho( 5 I 7 -> 5 S 2 ) 

P 91 = 2.59xlQ' 23 cm 3 /ps 





















8.2 SyssestiQiyJxtf-JEntHr^.Wark 


This thesis represents a comprehensive study of the spectroscopy and excitation 
dynamics of YLF:Tm,Ho at room temperature. A proper characterization of the 
temperature dependence of many of the phenomena discussed in this thesis has yet to be 
done. Phonons play an important role in the absorption, emission and energy transfer 
processes in Rare-Earth doped solid-state materials, and a better understanding of some of 
these processes would be gained by studies of the behavior at low temperatures, and at 
high temperature. This is especially true in the case of the Tm-Tm cross-relaxation and 
Tm->Ho transfer. It would be interesting to see if the models proposed in this thesis to 
describe the observed decay behavior remained valid at low temperatures and at high 
temperatures. 

Another area that provides a direction for future work is the use of more 
sofisticated techniques to probe into the energy transfer processes. Time resolved studies, 
for instance, can reveal information about the individual Stark levels participating in the 
transfer process by observing the fluorescence at various delay times after the excitation 
pulse. Time resolved studies of the Tm->Ho energy transfer would improve the 
understanding of the Tm<->Ho energy sharing. 

Concerning the models used to describe the experimental data in this thesis, there 
is always room for improvement. New experimental techniques lead to a greater finesse in 
the methodology that inevitably produces new physical insights into the processes being 
studied. The new physical insights allow for refinements in the present models or 
possibly entirely new models. 
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Appendix A Bfi cr-Lambert Law for a Uniaxial Crystal 


For a uniaxial crystal, such as UYF4, the polarization variation for the attenuation 
of an electromagnetic wave can be derived. We consider an incident electromagnetic wave 
propagating in an arbitrary direction described by the wave vector k. This is pictured in 
figure A-l. 



Figure A-l Coordinate System for Arbitrary direction of Propagation 

The angle between the wave vector k and the c-axis is 0 . The incoming E-vector 
is randomly oriented in the p-a plane, where the p-axis denotes the polarization axis. The 
orientation of E transverse to k is given by the angle xjr. For a beam (with a polarization 
described by y ) of intensity I 0 incident on a uniaxial crystal at some arbitrary angle 0 
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with respect to the c-axis, the transmitted light through the crystal can be written, after 
some geometrical considerations as: 

I(X;0,<p, y) = I Q jsin 2 y sin 2 0exp[-a Jt (A.K| + (cos 2 y + sin 2 y cos 2 0)exp[-a a (X)^]| 

(A-l) 

Note that this result is independent of the azimuthal angle <{> . For an unpolarized 
beam, we must average over the polarization angle y from 0 to 2 jc . The result is: 

I(X;0,y) = ±I 0 {sin 2 0exp[-a„(X)f] + (1+ cos 2 0)exp[-a CT (X)f]} (A-2) 


The spatially averaged transmitted intensity over the angles 0 ,<J> is: 


2 n it. 


| d<pj jsin 2 0exp[-a n (X)^] + (1 + cos 2 0)exp[-a a (X)f jjsin 0d0 
IW = ilo JL — 2 *r-; (A-3) 


2 it it 

J d<p J si 

0 0 


Performing this integration we find the result: 


IO-) = I 0 {^expl-a J (>.)fl+|exp[-a o (A.V]} (A-4) 


The absorption coefficients can be written in terms of the absorption cross 
sections Ojc=a K /N and o<j=<x<j/N. Thus, eq. 5-1 is obtained: 

I(^) = Io|^exp[-0 Jt (X,)N^] + ^exp[-o a (X)Nf]| (A-5) 
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Appendix B Least Squares Fitting Equation for Judd-Ofelt Parameters 


In order to find the set of Judd-Ofelt parameters that minimizes the sum of the 
squares of the differences between the theoretical and experimental linestrengths, begin 
with the matrix equation 5-5 for the theoretical linestrength: 


S = MQ 


0 - 1 ) 


This can be written in component form as: 


Sj = jMjiQi (B-2) 

i=l 


We want to minimize the expression: 

X 2 =t(Sj-Sf) (B-3) 

i=l 


where Sj E is the experimental linestrength. This is achieved by taking the derivative: 



_ d 


an 


k j=4 


XMjiQj- 


Li=l 



(B-4) 


and setting it equal to zero. This is easily done, and the result is: 


E M jki M jA = t MjtSf (B-5) 

j=l i=l H 
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The right hand side of eq. B-5 has the same form as eq. B-2 for the theoretical 

linestrength. The only difference is in the order of the indicies. So, recalling that the 
definition of the adjoint of a matrix, we can write = m£j . Eq. B-S becomes: 



(B-6) 


This substitution has the effect of turning the right hand side of eq. B-5 from a 
row matrix to a column matrix. In matrix form this is: 

M t MQ° = M*S e (B-7) 


Eq. 5-6 then follows from the above expression: 

£2° = (M t M)' 1 M t S E (B-8) 


which gives the set of Judd-Ofelt parameters Q° that minimizes eq. B-3. 


148 



Appendix C System Response Calibr ation for L uminescence Measurements 


It is possible to take into account the signal response of the experimental 
apparatus as a function of wavelength. Depending upon the spectrometer grating, 
detector, optics, polarizers, etc. the measured signal will have a wavelength dependent 
variation. For example, a grating has its highest efficiency at the Blaze wavelength, and 
becomes less efficient the further the wavelength is from the Blaze wavelength. Detectors 
have a characteristic responsivity that varies with wavelength. All of the effects that 
cause fluctuations in the signal level with changes in wavelength can be corrected for with 
the use of a lamp of known temperature. 

Suppose a luminescence spectrum is obtained. By introducing a correction 
function h(X.) we can write: 


Lmeas = (C-l) 

where L me as is the measured luminescence spectrum and L rea ] is the actual luminescence 
spectrum. The function h(X) takes care of the wavelength dependent changes in signal 
level due to the experimental apparatus. 

Similarly, we can obtain a spectrum of the temperature lamp. The temperature 
lamp simulates a blackbody source and will have a spectrum that is similar to a Planck 
distribution curve for the given temperature. Therefore we can write: 

^meas = h(X)P planet (C-2) 

Combining equations C-l and C-2 the real luminescence spectrum is found: 
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(C-3) 


^real 


planck - 


meas 


"meas 


In order to find the real luminescence spectrum, we must obtain the luminescence 
and lamp spectra under the same conditions and wavelength range. The Planck curve can 
be calculated if the temperature of the lamp filament is known. A convenient form for the 
power radiated by a blackbody per unit area per wavelength interval is: 


W X = 


Cl 

X 5 [exp(C 2 AT)-l] 


where, 

Cj = 2itc 2 h = 3.746x105 mW-cm 2 
C 2 = hc/k B = 1.438x107 nm-°K 
T = lamp temperature 
c = velocity of light 
h = Planck's constant 
kg = Boltzmann's constant 

All of the spectra required to calculate the real luminescence spectrum corrected 
for system response through eq. C-3 can be combined with a suitable computer program 
or a commercial graphing application such as Kaleidagraph or Cricket Graph. 
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Appendix D Ctosed Form Solution to Rate Equations for Tm->Ho Energy Transfer 


We begin with the linearized form (n 2 n 7 =0) of equations 7-14: 

~ an 2 + P n 7 (D-la) 

^l = -|2-+an 2 -Pn 7 (D-lb) 


We take as a trial solution: n 2 (t) = A 2 exp(-yt) and n 7 (t) = A 7 exp(-yt). The 
trial solutions are then substituted back into equations D-la and D-lb. The result is two 
equations: 


^_X- a ^A 2 + pA 7 = 0 (D-2a) 

(y-^-P)a 7 +“A 2 =0 (D-2b) 

Solving these equations for y leads to a quadratic equation: 

72 -(i + i +a+p Mvf7 + t + i) =0 (D ' 3) 


The solution to eq. D-3 is: 


^ = i(i + i + “ +l 


p) ± lJ(i + i +o+ p) 2 -fe + t + ^) v-v 
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This can be rearranged into the form: 


*=2(*r + ^ + “ + 



+ 2(a-(l)^4~ ^-j+(a+p) 2 (D-5) 


At this point we can make the reasonable assumption that the energy transfer 
rates are much faster than the lifetimes of the Tm,Ho manifolds. That is: 



(D-6) 


Equation D-5 can now be written as: 


= 1(4 


2^X 7 T 2 


+ 0C + 



(P-7) 


Invoking the binomial expansion (1+ x) 1/2 = 1+ nx + [n(n-l)/2]x 2 + ... and keeping 
only terms to first order: 


Y = 


h{i- 7 + ^- + 




a 


+ P + 



(D-8) 


This equation leads to two solutions: 



(D-9a) 


(D-9b) 
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The trial solutions are now modified to be a linear combination of exponentials for 


the two possible gammas: 


n 2 (t) = A 2 exp(-Yit) + B 2 exp(-y 2 t) (D-lOa) 

n 7 (t) = A 7 exp(-Yit) + B 7 exp(~Y 2 t) (D- 1 Ob) 

The exponential coefficients are found by substituting equations D-9 into 
equations D-10 and by considering the initial excitation conditions (@t=0 pumping Tm 
n 2 (0) - n 2 o , n 7 (0) =0). Making the assumption again that the energy transfer rates are 
much faster than the lifetimes of the manifolds and that the manifolds reach equilibrium 
quickly (x 2 = x 7 ) then the final result is: 

^ ' T) + (^Tp) ex P [ ~ ( “ + W (D - 11 a) 

= (^) exp( "‘ ,T) - (5Tp ) expHa + 1 m (D ' 1 1 b) 

which are the equations 7-14a and 7- 14b. 
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p 

We begin with writing the rate constants as a summation over the initial and final 
levels of Ho, and the initial and final levels of Tm for the Tm<— > Ho energy transfer: 

P28 =ZX w i*PiP* (E-la) 

jk it 

P71 = ZE w jkPjPk (E-lb) 

i l jk 

Writing the Boltzmann fractions out explicitly, these equations can written as: 

P28 ~ 7 . 7 0 ^ exp(~AEj / kgT)exp(— AE^ / kgT) (E-2a) 

28 it jk 

P71 = zt XX w jk exp(-AEj / k B T)exp(-AE k / kgT) (E- 2 b) 

7 1 i* jk 

The Stark level energies for an ion in the upper manifold can be written in terms of 
the zero line energy. Stark level energy of the lower manifold, and the energy difference 
between the upper and lower Stark levels. This procedure is similar to that discussed in 
chapter 2, section 6 . The result is: 

AEi = AEj + hvy - Eg 4 (E-3a) 

AE k = AE^ + hv k £ - Ezl (E-3b) 
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Inserting these in equations E-2 we find: 


P28 = -- (E ^ z ^ kBT) XX w i< ex p[(— A E j-A E /)/k B T]exp(— hv ij/k B T) (E-4a) 
28 U jk 

517 

P71 = e ~ (E |^ z / ^ T) w Jk cxp[(-AE j -AE ( )/ k B T]exp(-hv^ / k B T) (E-4b) 
71 U jk 

For resonant transitions we would expect Vy = Vyj, and w^ = wjy. So, the terms 

in summation in E-4a and E-4b are equal. Taking the ratio of E-4a and E-4b, then, we have 
Equation 7-16: 

(E-5) 

This equation is a general relation for any resonant sets of energy levels, and is 
applicable to other sets of energy levels, not just the ones discussed here. 
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